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Abstract 

We study a small quantum system (e.g. a simplified model for an atom or molecule) inter- 
acting with two bosonic or fermionic reservoirs (say, photon or phonon fields). We show that 
the combined system has a family of stationary states, parametrized by two numbers Ti, T2 
("reservoir temperatures"). If Ti 7^ T2, then these states are non-equilibrium, stationary states 
(NESS). In the latter case we show that they have nonvanishing heat fluxes and positive entropy 
production. Furthermore, we show that these states are dynamically asymptotically stable. The 
latter means that the evolution with an initial condition, normal with respect to any state where 
the reservoirs are in equilibria at temperatures Ti and T2, converges to the corresponding NESS. 
Our results are valid for the temperatures satisfying the bound min(ri,r2) > g^'^°' , where g 
is the coupling constant and < a < 1 is a power related to the infra-red behaviour of the 
coupling functions. 

1 Introduction 

The present paper is a contribution to rigorous quantum statistical mechanics. Key problems here 
are dynamical stability of equihbrium states, and characterization (if not the definition) and stabihty 
of non-equilibrium stationary states (NESS). 

While our understanding of the quantum equilibrium states, the subject of equilibrium statistical 
mechanics (see [HI HZl [THl ) , and the recent progress in proving their dynamical stabihty [51 
[TTl [TOl [22I [T4l [23] ) are satisfactory, results on non-equilibrium stationary states are just beginning 
to emerge. The problem is that we do not have a simple stationary characterization of NESS similar 
to the principle of maximum of entropy or the KMS characterization for equilibrium states. Thus 
it is remarkable that certain characterizations of NESS and their stability were recently shown 
for (idealized) particle systems coupled to Fermi reservoirs at high temperatures, min(Ti,r2) > 
C[lni]~^, in [21], for XY-chains [4] and for coupled Fermi reservoirs in [T2l|T5l|2l [1] . Here, Ti_2 are 
the temperatures of the reservoirs, and g is the coupling constant, which is assumed to be sufficiently 
small. 

There are two rigorous approaches to non-equilibrium, quantum statistical mechanics. One is 
based on scattering theory - wave (or M0ller) morphisms - and the other, on the theory of resonances 
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via complex deformations. In this paper we follow the second approach which we believe applies to 
a wider class of physical models. 

In this paper we establish a spectral characterization of the NESS and prove their dynamical 
stability for (idealized) particle systems coupled to two Bose reservoirs (e.g. photons or phonons) 
for reservoir temperatures satisfying min(ri,T2) > Cg^^", where C is a constant and < a < 1 
{a = jj^^yf, where fi > 1/2 is given in Condition (B) below). Our approach applies to an arbitrary 
finite number of bosonic or fermionic reservoirs; in the latter case, it gives an extension of the results 
of [21] to the temperature range mentioned above. Moreover, we develop a perturbation theory for 
NESS and use it to prove that the entropy production is strictly positive. 

An appropriate iteration of our estimates in the spirit of the spectral renormalization group 
of [3 O [71 [H] would give the above results for all temperatures. This extension will be presented 
elsewhere. 

Similarly to [21], we construct a NESS from a zero (non-degenerate) resonance eigenvector of a 
certain non-self-adjoint Liouville operator, K , acting on a positive temperature Hilbert space. The 
operator K is the analogue of a C-Liouvillean in the terminology for C* -dynamical systems 2li. 

To show dynamical stability of a NESS we have to establish certain long-time (ergodic) properties 
of the evolution, U{t), generated by K. The operator K does not belong to a class for which the 
evolution is a priori known to exist (e.g. a class of normal or accretive operators). To overcome 
this problem we establish a direct connection between desired ergodic properties of U{t) and certain 
spectral properties of a complex deformation. Kg, 9 G C^, oi K. For technical reasons we can use 
neither the complex deformations introduced in [18] nor those introduced in [8] but we combine both 
types, hence 9 is in rather than in C. (Such a combination was already mentioned in [S]). In 
order to establish the desired spectral characteristics of the operator family Kg , we use the method 
of the Feshbach map, as developed in [SKUlI]. 

The present paper suffers from the main weaknesses shared by all the works in the area, except, 
in some aspects, of [8]: 

(i) The particle system has a finite-dimensional state space; 

(ii) The restriction on the coupling fimctions is severe; 

(iii) Temperatures considered are high. 

To overcome the first limitation one would have to go beyond, or at least significantly extend, the 
present approach. The second limitation is due to use of translation analyticity (which in our case 
is combined with the dilation analyticity), see Remark 3 in Section 3. This analyticity is used in the 
present work in a single place - in controlling the nonsingular part of the resolvent of the operator 
Kg near the zero resonance pole by rendering this pole isolated and therefore the nonsingular part 
of the resolvent analytic (see the estimate (|8.8p ). Without it the zero resonance of the operator Kg 
is not separated from the continuous spectrum and sits exactly at a threshold of the latter. Hence 
to control the nonsingular part of the resolvent near the zero threshold becomes a delicate matter. 

The paper |8j has rather mild restrictions on the coupling functions due to using the dilation 
analyticity. Since [8] deals with the dynamics near equilibrium, the operator K in this case is 
self-adjoint and an analogue of Eqn (|8.8|) is obtained with help of an abstract spectral theory of 
self-adjoint operators. Furthermore, ^ handles arbitrary temperatures by employing the spectral 
renormalization group. In the present paper we take the first step in removing the high temperature 
restriction. To this end we use the Feshbach map of [7] . Already a single application of the Feshbach 
map considered in this paper improves the temperature bounds yielding the results mentioned above. 
We also set the stage for the iteration of this map - the spectral renormalization group method - 
which would remove the restriction on the temperature altogether. The iteration procedure will 
be carried out elsewhere. (Note that the works [TH [TH] deal with arbitrary temperatures, but the 
scattering approach they use seems to be inapplicable to the models considered in this paper.) 
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A more detailed outline of our approach and of the organization of the paper is given in Section 

2. 



2 Model and Approach 

We consider a system consisting of a particle system, described by a self-adjoint Hamiltonian Hp on 
a Hilbert space Hp, and two bosonic reservoirs, at inverse temperatures Pi and P2, described by the 
Hamiltonians Hri and Hr2 acting on Hilbert spaces Tiri and 7ir2, respectively. The full Hamiltonian 
is 

H:=Ho+gv, (2.1) 
acting on the tensor product space Ho Hp (E) Hri ® Hr2- Here 

Ho := Hpi»l(El + li» Hri<»l + l<»l<» Hr2 (2.2) 

is the unperturbed Hamiltonian, v is an operator on Ho describing the interaction and g e M is a 
coupling constant. 

For our key results we have to assume that Hp is finitc-dimcnsional, but some of our results hold 
for dimHp = 00. 

The operators Hrj describe free scalar (or vector, if wished) quantum fields on Hrj, the bosonic 
Fock spaces over the one-particle space L'^ {M."^ , k) , 

Hrj = j uj{k)a*{k)aj{k)d^k, (2.3) 

where a*{k) and aj{k) are creation and annihilation operators on Hrj and uj{k) — \k\ is the dispersion 
relation for relativistic massless bosons. The interaction operator is given by 

2 

v^^v.j with Vj = aj{Gj) + a*{Gj). (2.4) 

Its choice is motivated by standard models of particles interacting with the quantized electromagnetic 
field or with phonons. 

In (|2.4p . Gj : fc I— > Gj{k) is a map from R'^ into B{Hp), the algebra of bounded operators on Hp, 
and 

and a-(G,,:=«,(G,)-. (2.5) 

If the coupling operators Gj are such that 



^ y" (H- \k\-^) \\Gj{k)f dk is sufficiently small, (2.6) 

R3 



then the operator H is self-adjoint (see e.g. [5]). 

Now we set up a mathematical framework for non-equilibrium statistical mechanics. Operators 
on the Hilbert space Ho will be called observables. (Only certain self-adjoint operators on Ho 
are actually physical observables.) As an algebra of observables describing the system we take the 
C*-algebra 

A^B{Hp)(g>WiLl)(g>W{Ll), (2.7) 

where 21J(ig) denotes the Weyl CCR algebra over Ll := L'^{R^, (1 + \k\-^)d^k), i.e. the C*-algebra 
generated by the Weyl operators Wj{f) := e'*^^f\ (/),{/) := ^ (a*(/) + a^if)), with / G L^, see 
e.g. [H]. States of the system are positive linear ('expectation') functional ip on the algebra A, 
normalized as V'(l) — 1- 
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The reason we chose A rather than B{'Ho) is that the algebra A supports states in which each 
reservoir is at a thermal equilibrium at its own temperature. More precisely, consider the evolution 
for the j-th reservoir given by 

alj{A) ■.^eJ"'-^'Ae-'"'-^K (2.8) 

Then there are stationary states on the j-th reservoir algebra of observables, W{Lq), which describe 
(single-phase) thermal equilibria. These states are parametrized by the inverse temperature (3j = 
1 /Tj and their generating functional is given by 

4f (W^.(/))=exp|-l |^^^J^|/(fc)Pd3,|. (2.9) 

The choice of the space Lq above is dictated by the need to have the r.h.s. of this functional finite. 
These states are characterized by the KMS condition and are called the (a*^ , /3j)-KMS states. 

Remark. It is convenient to define states ip on products a'^ifi) ■ ■ ■a^{fn) of the creation and 
annihilation operators, where denotes either a or a* . This is done using derivatives ds^. of its 
values on the Weyl operators W[sifi) . . . W{snfn) (see 0, Section 5.2.3 and (I2.15P ). 

Consider states (on A) of the form 

uo:=ujp®uji{'^ (2.10) 

where Wp is a state of the particle system and w^f^ is the (a* j, /3)-KMS state of the z-th reservoir. 
The set of states which are normal w.r.t. lhq is the same for any choice of LOp. A state normal w.r.t. 
Wo will be called a /3i/32-normal state. 

In the particular case Wp(-) = Tr(e~'^p'^'' •)/Tr(e~'^''^'') we call wq a reference state. 

The Hamiltonian H generates the dynamics of observables A G B{'Ho) according to the rule 

A ^ a\A) e'^'Ae''"' . (2.11) 

Eqn (|2.1ip defines a group of *-automorphisms of B{TLq). However, a' does not map the subalgebra 
A C B{Ti,Q) into itself, so (|2.1ip does not define a dynamics on A. To circumvent this problem 
we define the interacting evolution of a class of states on A by using the Araki-Dyson expansion. 
Namely, we define the evolution of a state ip on A which is normal w.r.t. luq by 

^\A):^ lim y^iigr / dt, ■ ■ ■ / V"^' ■*™(A), (2.12) 

where the term with m = is ^(ao(v4)), and, for m > 1, 

^Mi,...,t™ (^) ,^ ^ ([^t™ (^^)^ . . . [^U (^)] ...]). 

Here, = w* G ,4 is an approximating sequence for the operator v, satisfying the relation 

lim c^o(^*« - v*)ivr, - v)A) = 0, (2.13) 

n — >oo 

for all A polynomials in a*{f), j = 1, 2, / G Lq. Such a sequence is constructed as follows. Let {e„i} 
be an orthonormal basis in Lq. We define the approximate creation operators 

= E Gs)blx{em), (2.14) 

m— 1 

where n = (A,^), and, for any / G L'^(M'^) and A > 0, 

b*,xif) ■■= ^ {Wjif/X) - 1 - ^W,{^f/X) + zl} . (2.15) 
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Similarly we define the approximate annihilation operators aj^n{Gj). These operators belong to A. 
Via the above construction we obtain the family of interactions Vn which belongs to A and, as can 
be easily shown using l|2.9p . satisfies (|2.13p . 

We show in Appendix [X] that under condition (|2.13p the integrands on the r.h.s. of (|2.12p are 
continuous functions in ii, . . . that the series is absolutely convergent and that the limit exists 
and is independent of the approximating sequence Vn- 

Our goal is to understand stationary states of the interacting system originating from [3i[32- 
normal states either by a perturbation theory or through an ergodic limit of the full evolution a*. 
These states are not equilibrium (KMS) states states. They will be called non- equilibrium stationary 
states or NESS for short. Their main feature is that the energy (heat) fluxes between the reservoirs 
and the particle system do not vanish. 

Assuming certain smoothness and smallness conditions on the coupling operators gGj{k) and 
assuming that the particle system is effectively coupled to the reservoirs, we show that, starting 
initially in any /3i/32-normal state ■0, the system converges, under the evolution a*, to a state 77: 

■0* — >ri ast^oo. (2.16) 

The convergence (|2.16p is understood in the weak* sense on the sub-C* -algebra of "analytic observ- 
ables" 

Ai = B{Hp) ^ WiD''''^^) (g) 211(1?''"^'). (2.17) 

Here, *2II(£>'*"^') is the Weyl CCR algebra over the dense set l?'*"^' c which we define in Ap- 
pendix [Cl Roughly speaking, D^-nai consists of vectors from the space nb>o ^"''''^'-^o which have 
some analyticity properties in |fc| and a certain behaviour at fc = 0. The density of D^'^^^ C Lq 
implies that Ai is strongly dense in A. The construction of the state 77 and the proof of its stability, 
(|2.16p . will rely on the theory of resonances for the evolution ^ ^ 

As mentioned in the introduction, so far, we do not have a simple characterization of NESS. How- 
ever, there is a key physical quantity which differentiates between equilibrium and non-equilibrium 
stationary states - the collection of heat fluxes. In our case, the heat flux, or more precisely the 
heat flow rate (i.e. the energy flow rate due to thermal contact), j =^ 1, 2, into the j-th reservoir 
is given by 

^r-= §i\t=o^\Hr,), (2.18) 
and the heat flux, ipQ, into the particle system is defined as 

<^o:=|L=o"*(^p)- (2-19) 

The heat fluxes can be combined into a single quantity - the entropy production. Motivated 
by the second law of thermodynamics {dS = '^(ijdQj) we introduce the observable of entropy 
production (rate) as 

2 

s:=^/3,0,, (2.20) 

j=o 

where, for notational convenience, we write f3o := f3p. The entropy production, EP(uj), in a state co 
is defined as (see [36l Ell EHl |32l [33l [Ml ESI [13 [20l [2^ 

EP{uj) ^ uj{s) . (2.21) 

Since s is not a bounded operator, we have to use an approximation procedure similar to the one 
mentioned in the remark after (j2.9p in order to define the r.h.s. of (j2.2ip for sufficiently regular 
states. 
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The entropy production EP{rj) of the NESS is independent of the particle state ujp entering 
Definition (|2.10|) of the state of the decoupled system, since 77 is independent of ujp. Notice that 
vii'o) = dt\^_^ii{a* {Hp)^ , since Hp is a proper observable (here we assume Hp e B{Hp)) and 77 is a 
continuous and stationary state. Hence 

r/(0o) = . (2.22) 

2 2 

Therefore, writing v{(l>j) = 'ni'l^j) = ?7(^t|t=oa*(-ffo)) = -dt\t^QViaH'")) = 0, we obtain 
3=1 3=0 

2 

^77(0,) -0. (2.23) 

Observe that the zero total flow relation (|2.23[) and Definition (|2.20p for the entropy production 
rate imply that 

i?P(77) = (/3i-/32)77(0i) . (2.24) 
Thus, the relation EP{r]) > is equivalent to 

77(01 ) > whenever Ta > Ti , (2.25) 

where Tj = [3~^ is the temperature of the j-th reservoir. In other words, in the state 77 the energy 
flows from the hotter to the colder reservoir. 

A general resuh due to ^ shows that EP{uj) > for any NESS uj. We show that for the NESS 

EP{7^) > iff Pi ^132, 

see Theorem 13.21 and Section [12] for a precise statement of this result. Moreover, we develop a 
perturbation theory for the NESS and compute EPijf) in leading order in the coupling constant g. 

Let us outhne the main steps of our proof of the convergence ()2.16|) (c.f. [21]). We pass 
to the Araki- Woods GNS representation of (^, wq), with wq of the form (|2.10p and ojp{A) := 

{A,LOo) {H,7T,no), 

where Ti., n and are a Hilbert space, a representation of the algebra A by bounded operators on 
7i, and a cyclic element in Ti. (meaning that Tr{A)flQ = Ti) s.t. 

LOoiA) = (Oo,7r(A)r!o). 

The GNS representation provides us with a Hilbert space framework which we use to convert the 
dynamical problem described above into a spectral problem for a certain non-self-adjoint operator 
K on the Hilbert space H. With the free evolution ao(A) :— ^i-tHo j^^-itHo Q^^g associates the unitary 
one-parameter group Uo{t) — e**^" on Ti, s.t. 

7r(a* (A)) = Uo{t)n{A)Uo{t)-^ (2.26) 

and Ua{t)il,o = fio- Define the operator L^^^ := La + g'!T{v) on the dense domain Dom(Lo) H 
Dom(7r(w)). Here tt{v) can be defined either using explicit formulae for tt in the Araki-Woods 
representation given below or by using the approximation, Vn G A, for the operator v constructed 
above. By the Glimm-Jaffe-Nelson commutator theorem the operator L^^^ is essentially self-adjoint; 
we denote its self-adjoint closure again by the same symbol L^^\ The operator L^^^ induces the 
one-parameter group cr* on tt{A)" , the weak closure of tt{A), 

f7*(B) :=e'*^*"Be-**^'" (2.27) 



6 



for any B S tt{A)" . Let ■0 be a state on the algebra A normal w.r.t. ojq, i.e. 

i:{A) = Tr(p7r(A)) (2.28) 

for some positive trace class operator p on 7i of trace one. It is shown in Appendix [K\ that for V' as 
above the limit on the r.h.s. of (j2.12|) exists and equals 

V'*(A)=Tr(pa*(^(A))). (2.29) 

In particular, the limit is independent of the choice of the approximating family w„. 

Due to (|2.29p the dynamics on normal states, defined in (|2.12p . gives rise to the dynamics on the 
Hilbert space 7i, determined by a one-parameter group U{t), satisfying 

Tr{pU(t)TT(A)U{t)-^) ^i/{A), VAeA. (2.30) 

Due to the fact that the von Neumann algebra tt{A)" has a large commutant (which isomorphic 
to Tr{A)" , as is known from Tomita-Takesaki theory), relation ()2.30p does not define U{t) uniquely; 
however, if we impose in addition to (|2.30|1 the invariance condition 

Uit)n = fl, (2.31) 

where 51 is a fixed cyclic and separating vector, then U{t) is uniquely determined. (The vector f2 
is called cyclic if 7T{A)il. is dense in H and separating if 7r(^)'Sl is dense in H, the prime denoting 
the commuant.) If H. were the vector representing an equilibrium state then U{t) satisfying (|2.30p 
and (|2.3ip would be a unitary group. In the non-equilibrium case /3i (32, one can see that (|2.3ip 
cannot be satisfied for a unitary U(t) implementing the dynamics as in ((230)) . For technical reasons, 
we choose U{t) to satisfy p.3ip for a convenient vector fi, rather than to be unitary (c.f. [21]). 

We will show that U(t) is strongly differentiable on a dense set of vectors and we will calculate 
explicitly its generator, K :— —i-g^U{t)\t=Q. In the non-equilibrium situation K* ^ K {U{t) is not 
unitary!) and (|2.3ip implies that Kfl — 0. The main effort of our analysis is to derive enough 
spectral information on the operator K to enable us to show (|2.16p and to identify the NESS with 

ri{A) ^ {n*,TT{A)n) , (2.32) 

where fi* is a zero resonance of the operator K*: K*n* = (in the sense of distributions) and 
52* S ^'anai' ^'^^ appropriate dense set I?anai C Ti., and A are such that TT{A)il G V^nai- 

In order to obtain rather subtle spectral information on the operator K, and to give a precise 
meaning to expression (|2.32p . we develop a new type of spectral deformation, K l—^ Kg, with a 
spectral deformation parameter G C^, in combination with an application of a Feshbach map 
acting on Kg. 

In conclusion of this outline we present here the GNS triple provided by the Araki- Woods con- 
struction, which forms a mathematical framework for our analysis (see [8j [181 IS] for details and 
O [II] for original papers). In the Araki- Woods GNS representation the (positive temperature) 
Hilbert space is given by 

n^np®n\ (2.33) 

where W ^Hp® Hp and TV = TT^ ® with 

rC^ ='Hr3®nrj. (2.34) 

We denote by afj{f) (resp., afj{f)) the creation and annihilation operators which act on the left 
(resp., right) factor of (|2.34p . They are related to the zero temperature creation and annihilation 
operators a'^if) by 

7r(a,(/)) = a,,(yrT^ /) + al^{^J) (2.35) 
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and 

= %,(Vp7/) + ^rji^l+Jj f) (2.36) 

where pj = Pj{k) = (e^j^C^') — 1)^^ with Ljj{k) = |fc|. Finally, we denote Vl^ :— flri <8) 5^r2, where 
ilrj ■— ^rj,e ® ^rj.r are the vacua in . Thus, fir is the vacuum in Ti,^ . 

Definition (|2.10p and our choice of ujp made at the beginning of this section imply that 



[^.e-'3p-Ej]i/2 



no = np(g) nr with np = nll'-'> = ^^f^^ J,,,^ \ (2-37) 



where Ej and ^pj are the eigenvalues and normalized eigenvectors of Hp. 

The self-adjoint operator Lq generating the free evolution ?7o(<) defined in (|2.26p is of the form 
Lo = Lp ® V + I'P ® Lr with Lj. = ^rj- The operator Lp has the standard form 



Lp ^ Hpig lp ~ lp ig Hp, 

and 

Lrj = J uj{k) [af j{k)aej{k) — a* j{k)arj{k)) Sk. 

The operator K can be written as i^T = Lo + g{y — W) with = tt{v) and = ti'{w) with w a 
non-self-adjoint operator obtained by a simple transformation of v. 

A standard argument shows that the spectrum of the operator Lq fills the axis M with the 
thresholds and eigenvalues located at u{Lp) and with an eigenvalue of multiplicity at least dimiJp 
and at most (dim i/p)^ (depending on the degeneracy of the spectrum of Lp). A priori we do not 
know anything about the spectrum of the non-self-adjoint operator K besides the fact that it has an 
eigenvalue 0. For all we know its spectrum might fill in the entire complex plane! Thus understanding 
the evolution generated by the operator ii' is a subtle matter. 

This paper is organized as follows. In Section 3 we give a precise formulation of our assump- 
tions, state the results and discuss assumptions and results. In Section 4 we present the Hilbert 
space framework and we define the vector f2 and the evolution U{t) and in Section 5 we describe 
the generator K. In Section 6 we introduce the complex deformation Kq of K and establish the 
connection between the resolvents of K and Kq. In Section 7 we establish the spectral properties of 
Kq which we then use in Section 8 to express the dynamics in terms of an integral over the resolvent 
of Kq. In Section 9 we prove our first main result, the existence and explicit form of the NESS, 
and its dynamical stability. Section 10 contains a supplementary result on the structure of the level 
shift operators, used in Section 7. In Section 11 we develop a perturbation theory for NESS and in 
Section 12 we prove the positivity of the entropy production. Finally, in Appendices lAl lP] we collect 
some technical results. 



3 Assumptions and Results 

In order to state assumption (B) below, it is practical to define the map 7 : L^(R'^) L^(R x S"^), 

( f\( \ /n / /('"^)' w > 0, , . 

(7/)(u, cr = V \ -Fr \ r. (3-1) 

Let jQ{u) = e'^"sn(«)y _^ ^ fo^. q ^ (5^7-) e and w e R (see jM!])) and define (7e/)(u,cr) = 
(7/)(je(w), cr), for / e i^(R x 5^), 6* G R^. The maps 7 and 79 have obvious extensions to operator 
valued functions. 

(A) Ultraviolet cut-off. J \\Gjik)fe''\''\^[l + \k\-^](fk < 00 for some a > 0. 
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(B) Analyticity. For j = 1,2 and every fixed {u, tr) S R x 5^, the maps 

9^heG,){u,a) (3.2) 
from to the bounded operators on Tip have analytic continuations to 

{{d,T)eC^\\ImS\<d„,\T\<To} , (3.3) 

for some Sq, tq > 0, < where /3 — max{(3i , P2) ■ Moreover, 

1/2 

/ \/\u\ + 1 ^fl.U,l/9^ \ , 

79 



i/=l/2,/j 





2 


j (u,ct) 


duda 



< 00, (3.4) 



for some fixed ^ > 1/2 and where (5/3^ = /3 — For future references, := {Sq, tq). 

(C) Non-Degeneracy of the Particle System. We have dim Tip = N < 00 and the Hamiltonian Hp 
has non-degenerate spectrum {En}^^Q- 

(D) Fermi Golden Rule Condition. We have, for j = 1,2, 



70 := min / (5(|fc| - S„„)|G,(fc)„„|^d^fc > 0, (3.5) 

0<n<m<A^— 1 J 

where E'mn = Em — En, Gj{k)mn '■— (Vmj Gj{k)ipn), the (p„ are normalized eigenvectors of Hp 
corresponding to the eigenvalues E„, and S is the Dirac delta distribution. 

(E) Either dim Tip = 2, or dimTYp > 3 and the inverse temperatures satisfy 

|/3i - /?2| < c, or |/3i - /32| > C and min^ (3^ > C , (3.6) 

where < c, C < 00 are constants depending only on the interaction Gi^2- 

Remarks. 1) The map (|3.ip has the following origin. In the positive-temperature representation 
of the CCR (the Araki- Woods representation on a suitable Hilbert space, see Appendix A), the 
interaction term Vj is represented by aj {jp^ Gj ) -I- a* {jp. Gj ) , where 



2) The Condition (A) can be removed at expense of a slightly more involved proof of existence 
of the operator r(z) in Lemma [4.21 see the remark after the proof of Lemma [4.21 However, since 
the ultraviolet behavior of the coupling functions is inessential in the models considered, we choose 
a stronger assumption over extra steps in the corresponding proof. 

3) A class of interactions satisfying Conditions (A) and (B) is given by Gj{k) = g{\k\)G, where 
g{u) = u"e-"', with m > 0, a = n -I- 1/2, n 0, 1, 2, . . ., and G = G* G BiHp). A straightforward 
estimate gives that the norms p.4p have the bound 

l|G,ll,,^,,,,<c(l + e(^ft)^/4)||G||, (3.8) 

provided < a + 1, where the constant G does not depend on the inverse temperatures, nor on 
varying in any compact set (compare this with the bound (4.13) of [21 ). The restriction a = n-f- 1/2 
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with n = 0, 1, 2, . . . comes from the requirement of translation analyticity (the r-component of 9), 
which appears also in [21\ . 

4) The condition tq/ cos(5o <2-k/ (3 after (|3.3p guarantees that the square root in (13. 7p is analytic 
in translations m m + r. 

5) What we need in our analysis is that the level shift operator Aq, the N x N matrix defined in 
(j7.ip , has a spectral gap at zero which is bounded below by a strictly positive constant independent 
of the temperatures. Condition (E) ensures this property. If one can show the desired property of 
the gap by other means then Condition (E) can be dropped. 

Let 

a := min{|A -^1 \ X, fi e a{Hp),X^ fj.} (3.9) 

and define 

-1/2 ^ «-l/2^ 



(l+/3l + (^2 )niax sup II II 5^3^, 1/2,, 
J \9\<do 



go :=Ca'/^ sin(5o) 

where C is a constant depending only on tan(5o, and SPj ^ f3 — Pj, SPp — |/3 — /3p|, and set 

gi := min ((.go)'/", [rmn{/3^\ P^')]^) , (3-10) 

where a = and /i > 1/2 is given in (|3.4p . 

The main results of this paper are given in the following theorems, where by a "state" on 
a subalgebra (which is not necessarily a C*-subalgebra), we mean a positive normalized linear 
functional. 

Theorem 3.1 Assume conditions (A) ~ (E) are obeyed for some < /?i,/?2 < oo, fi > 1, and let 
(3 = max(^i,/32). 

< |.g| < gi then there is a stationary state rj — rjp-^p^, defined on a strongly dense subalgebra 
Ai of A (see {2. 7^ ), satisfying 

ip^-^rj, t^oo (3.11) 

for any (3if32-normal initial state ip. rj is continuous in the norm of A. The convergence is in the 
weak* sense (i.e., pointwise for each A g Ai). For A G Ai, r](A) is analytic in g. 

Remark. 6) Our analysis shows that the NESS is actually defined on a bigger (but somewhat less 
explicit) Banach space of operators 2 Ai (see (|9.1ip ). and the convergence to the NESS, (|3.1ip . 
holds on On Ai one can introduce a "deformation norm" ||| • ||| < || • ||, see (|9.10p . such that 
in this norm, the convergence in p. lip is uniform, sup^g_4^ |V'*(^) ^ ^(^)l/lll^lll ~^ 0- Moreover, 
on Ai , the convergence is exponentially fast for initial conditions t/j in a dense set (in the topology 
of bounded linear functionals on A) - this set is the convex hull of vector states with deformation 
analytic vectors. 

Theorem 3.2 Assume that the conditions of Theorem \3.1\ are satisfied and let r]f3-^f}^ and g be as in 
Theorem \S.l\ Let /3i ^ (32, and let g and \Pi — P2\ be sufficiently small (independently). If either 
Gi = G2 or dimT^p = 2, then EP{r]p^[}.^) > 0. 

Our analysis gives a stronger result than the one presented in Theorem l3.2l Namely, for > 3/2, 
we show that EP{rip^f3^) > 0, provided o{g°)0{d/3) < i , where 5(3 = |/3i - ;92| (see Theorem nTT]) . 
Here, rj' depends on the inverse temperatures and the coupling functions and is given by 



where Eji = Ej—Ei, gji{E)'^ — J^^ d^k\ {ipi, Gi{k)(pj) \'^5{Eji — uj) (see Condition (C)). The numbers 
7-,- > are the coordinates (in the basis {ipj ® ipj} of Null(Lp)) of the unique vector Q* in the kernel 



10 



of the adjoint level shift operator Ag, at the value f3p = (and normalized as J2j Ij — V^)- (The 
operator Aq is defined in Section [TOl) 

By general arguments one can prove that rj' > for /3i > /32, < for (3i < 02, and 77' = if 
/3i = ^2- We also show that 77' > for /3i > /?2 and 77' < for f3i < P2, for all Pi, (32, except possibly 
for finitely many values in any compact set, see the remark at the end of Section [T2l 

The dependence of rj' on Sj3 is determined by the coordinates 7^ . We compute those in the cases 
when Gi = G2 and dim Tip = 2 (see the proof of Theorem ll2.H and equation (|11.22p . respectively). 

Remarks. 7) Using Araki's theory of perturbation of KMS states, one shows that if the temper- 
atures of both reservoirs are equal then the limit state is an equilibrium state and has zero entropy 
production. Non-existence of equilibrium states for (3i ^ /?2 has been shown in [25] . 

8) For a model with fermionic reservoirs, using a sufficiently fast convergence rate in (|3.11|) (e.g. 
0[ir'^) with a > 1 suffices) and the fact that r] is not a normal state for /3i 7^ /3i, it has been shown 
by an abstract argument that EP{r\p^ji^ > 0, provided — /?2| > Cg for some C > (see [H]). 
Instead of this indirect derivation we compute EP{rii}-^[}^) to the leading order in g and derive the 
results stated in the theorem. 

9) The condition Gi = G2 can be relaxed to Gi — G2 being small in a suitable sense. 



4 Spectral Theory of NESS 

In this section we outline a spectral theory of NESS applicable to Bose and Fermi reservoirs. Our 
approach follows the one developed for the Fermi reservoirs in [21]. Fix a state, ujq, of the form 
([2l0]) with ujp{A) Ti{e-'^''"''A)/TT{e-f^''"''). 

In this and the next section we use Condition (A) and 

{Hp + i)^^ is of a trace class. (4-1) 

In particular, the 0-analyticity of the coupling functions and the finitcness of the dimension of the 
particle space are not required. 

We pass to the Araki- Woods GNS representation for the unperturbed system: 

{A,LUo) ^ (H,7r,r!g) 

where Ti, tt and Og are a Hilbert space, a representation of the algebra A by bounded operators on 
H, and a cyclic element in Ti. s.t. 

c^o(^) = {no,7TiA)no)- 

There is also an anti-linear representation, tt', of the algebra A in bounded operators on the space 
n, s.t. tt' commutes with tt (i.e. [tt' {A),tt{B)] = V^, B S A), and tt'{A)Qq = H. 
The full dynamics is implemented by a one-parameter group U(t) satisfying 

U(t)Bn = a\B)n, VB e TTiA), (4.2) 

where is a cyclic and separating vector for tt(A) to be specified below and where 

a*(i?) :=e"^'''i?e-'*^''' (4.3) 

for any B £ tt{A) with, recall, L^^-' := Lq + gTT{v) and Lq is defined on the line before p.26p (see 
also the paragraph after (I2.37|) ). Observe that ()4.2|) implies that 

U{t)VL = n. (4.4) 

If the state oj corresponding to f2 is stationary (i.e. lj* = uf) then U{t) comes out to be unitary. In 
our situation we expect that there is no wg-normal stationary state and U{t) will be a non- unitary 
group. 
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We pick the vector as follows. Let (3 = niaxj=i^2 f3j- We define 

r!:-e-^^'"/2f^„/||e-/3i<^'/2f^^||. (4.5) 

The facts that O is well defined, i.e., that Qq G Dom(e"''-^' '/^), and that O is cyclic and separating, 
are established in Proposition 14.11 at the end of this section. 

The family U{t) is not unitary since cu :~ (ri,7r( • )Q) is not stationary: 

{U{t)Tr{A)n,U{t)Tr{B)n) = {n,a*{Tr{A* B))n) ^ uj\A* B) 

^ uj{A*B) ^ {TT{A)n,TT{B)n) , (4.6) 

for some A,B,t. Let now ip be an WQ-normal state corresponding to the vector Q^l E H (i.e. 
i^iA) = {Qn,TT{A)Qn)), where Q € n'{A) then 

4>\A) = {Qn,a\7r{A))Qn) . (4.7) 

Due to Eqns (jA.4ll and crQ(7r(A)) :— tt (0:0(7!)), and due to the convergence of (T*^^(7r(A)) to iT*(7r(y4.)) 
established after (jA.4p . we see that the operator Q commutes with (T*(7r(A)). Using this together 
with and ((i^ we arrive at 

^b\A) ^ {Q*Qn,U{t)7:{A)n) . (4.8) 

This key formula, due to [H] , connects the long time behaviour of {A) with spectral properties of 
U{t) or its generator. We explain what this means. 

Assume we can show that, for a certain class of cj) and and as t ^ 00, 

{(t>,Um)^{4>,P^) (4.9) 

where P is the eigenprojection on the fixed point subspace of U{t) (i.e. U{t)P — PU{t) — P), which 
we assume for a moment to exist. Relations (14.81) and (I4.9P imply 

ip^^ri, t-yoo, (4.10) 

where the state rj is defined (on an appropriate set of observables) by 

T]{A) -.^ (Q*Qn,Pn{A)n) . (4.11) 



We will show below that (|4.9p holds for some unbounded projection operator P. To understand the 
structure of this operator, we proceed as follows. 

We will show that U{t) is strongly differentiable on a dense set and we will compute its non-self- 
adjoint generator, K := ~i-g^\t=oU(t), which satisfies (see (14.41) ) 

KQ = 0. (4.12) 

The operator P is the eigenprojection onto the eigenspace of K associated with the eigenvalue 
(i.e. KP = PK = 0). We show that dimP = 1 and 

p = \n){n*\ (4.13) 

for some fi* ^ H, satisfying K*n* = in a weak sense (fi* G ^'anai' where Panai = ^im0>o^om{Ug) 
with the family Ug defined in Section [5]). Understanding the nature of the vector Q* , which we call 
the NESS vector, is a goal of our analysis. 

Substituting (|4.13p into (|4.1ip and using that {Q*Qn, n) = WQnf = V'(i) = i, we obtain 

77(A) = 71(^)17). (4.14) 



12 



Since ^* ^ Ti. the state rj is not normal but it is well defined for a dense set of observables. The 
question now is what is fl*7 The answer, provided in subsequent sections, is that f2* is a resonance 
of K*. 

In the following sections we construct a mathematical framework which provides meaningful 
expressions replacing formal ones, (|4.9p - (|4.14p . and with the help of which we can prove the conver- 
gence (|4.10|) . 

Remark. Evolution groups and their generators given by conditions of the type of (|4.4p (or 
(|4.12|) ) were introduced in [21], where the group U{t) is specified by the condition U{t)flQ = ^Iq. 
where flo is the unperturbed vector ("vacuum") introduced in (|2.37p above. However, an analysis 
of the operator K (see Sections |6] and [71) defined this way requires, instead of Condition (B), the 
condition obtained from Condition (B) by replacing the weig 
an additional restriction on the temperatures of the form 

g<cmin{e"*}, i.e., minT^ > c[ln(l/5)]~\ (4.15) 

3 j 

Using in i2lj the vector instead of f2o would improve this bound to miuj Tj > eg. 

Now we proceed to the main technical result of this section - the proof of the existence of the 
vector D, and establishing its properties mentioned and used above. 

Proposition 4.1 e Dom(e-'^^*"/2) and the veetor fl := e-'3^'"/2r2o/||e"'^^*"/^floll is cyclic 
and separating for the von Neumann algebra 7r(^)". 

We begin with with some preliminary technical results. To do manipulations with unbounded 
operators we use the dense subset, T>, of our Hilbert space, H, defined by 

V 7r{A)no, (4.16) 

where A — B{T-Cp) ^Vi ^ 7^2- Here Vj is the polynomial algebra generated by the annihilation and 
creation operators, aj{f) and a*{f ), of the j-th reservoir acting on Tirj with / e L^(R'^,e° '''"I (1 + 
\k\~^)d'^k) for some a' > 0. I? is a subset of 

J^o := {i^ eH \ Pnip = ip, for some n < oo} , (4.17) 

where P„ is the spectral projection of the self-adjoint operator N = max{A''i, A^2} '■— 
/ max(A, ^)di?Ari (A) (g) dEM^iilA associated to the interval [0, n]. Here, Nj :— 

J a*i j{k)aej{k) + a* ,j{k)arj{k) d^k is the number operator of reservoir j (see also (|6.3p ). 

J-'o is commonly called the finite-particle subspaee. 

Let A and J be the modular operator and modular conjugation associated with the pair 
{^{A)" , VIq} and let be defined by 

= ® all' ® "^2*- (4-18) 
The vector ilo defines a (k, 1)-KMS state and one can show that 

A-"AA** =7r(K*(A)) VAeA (4.19) 

Using this one computes that A — e^^, where L = /3pLp + fir\i^r\ + Pr2Lr2, and, in particular. A** 
commutes with e^*^° . Observe that 

grLop ^ 2? and A^V = V Vt e M. (4.20) 

In the sequel, an important role is played by the operators 

T{z) e-^^'"/2g.Lo/2 (4 21) 

(defined as products of two, in general, unbonded operators). We study properties of these operators 
in the following lemma. 
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Lemma 4.2 The following statements hold for all z £ <C and all ip, ip € T): 

V C Dom(r(z)), (4.22) 

z V{z)lp is entire, (4.23) 
^ r(z)e-^-^«/2 o„ (4 24) 

(Br(z)V', T{-z)ip) = {B^l: ,Lp) VB e ti{A)' , (4.25) 

{JAT{z)na,if) - (A*f7o , r(z) A^/^) , VA e 7r(yl)", (4.26) 

T{z)V is dense. (4.27) 

Remark. 'k{A)' is the von Neumann algebra ir' {A)" , the weak closure of it' [A). 
Proof of Proposition Let T = r(/3), where r(z) is defined in (g^. Since e"'^«17o = l^o, the 

property fio G Doni(e^^^' '/^) is equivalent to fio G Dom(r) (see (|4.24l) ). which is proven in Lemma 
14.21 relation (|4.22p . Hence the vector Vl exists. 

For the cyclicity it is enough to show that if i? G i^{A)' and BQ, = then B ~ Q. Let 17' := 
g-/3LW/2nQ ^ rf^o- By Km . y^p € V, = {Br{^)no,^i-P)y^) = (^^^o This implies 

-BOo = 0. Since JIq is separating for 7r(y^)' we have that B = 0. 

Now we show that ft is separating for tt{A)" . Let A e T^iA)" be such that AfJ = 0. The relation 
Arrio — and equation (|4.26|) imply that 

= {JATflo, if) = I^A*VLa , FA^/V) , (4-28) 

for any ip <£ T> C J-q. Now A^^^V — V and TV is dense (as is shown in (|4.27[) '). so we have that 
FA^/^p is dense and it follows from (|4:28l) that A*no = 0. Since Qq is separating this implies that 
A — Q and therefore is separating. ■ 



Proof of Lemma \4.S\ We first show that for all ip £ J-q the following formal Dyson expansion of 
the operator V{z) is well defined: 



m>0 



f{z)p:= ^ (^) fdr,--- dr„y,^^-'\V)---ar'\V)^. (4.29) 



The integrals on the r.h.s. are understood as strong limits of Riemann sums. Due to the UV-cut-off 
Condition (A), the transformation a™{V) — e^'^^"Ve'^^° is well defined and strongly analytic on 
To for w £ C. Since (p E J-q, there is a i/q s.t. Piy„(p = p (see also (14.171) ). Because each interaction 
V can increase the number of particles in each reservoir by at most one, we can write the integrand 
of (j4:29l) as 

al^--/'{V)P,,+r.-i a^^^'"-/'(y)P.„+™-2 • • • al^^'/\V)P,,^. 

Since cr™(V^)Pfc are bounded operators, the integrand on the r.h.s. of (|4.29D belong to H. Moreover, 
it is strongly continuous in ri, . . . , r^. 

Let us first show that the series in (I4.29P converges absolutely, for all values of 5, z e C. We 
use the bound ||iTQ"'(l^)Piy|| < C(lmw){i' + 1)^/^ which follows in a standard way from the explicit 
expression of V. The constant is given by 
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It follows that the norm of the m— th term in the series of ()4.29p has the upper bound 

^^^^^'^li'^/^''^ — \f^^^^^' Consequently, the series converges absolutely, for any z e C, and for 
any ^ Tq. 

Next, we show that T) C Dom(r(2:)), and that T{z)tf — V{z)if for all (p E V. It suffices to 
establish that for any ip E T), 

for aU V' e H s.t. V = /(i^^^)V' for some / € C^(]R) (such form a core for e /^). Indeed, this 
would show that e^^^l'^ip S Dom(e~^^'*''/^) and therefore Lp £ Dom(r(z)). Equation (|4.30p can be 
shown e.g. using the analyticity of both sides in z, and the fact that the equation holds for z e iW. 
Indeed, in the latter case r(z) are bounded operators and the Dyson series expansion (|4.29p is valid 
for them. In particular, 2? C Dom(r(2:)). Moving e~^^' '/^ in (|4.30p to the right factor proves that 
T{z)ip = T{z)ip, for all (p £ V. This also shows (|4.23p . since the series (|4.29p is entire in z. Thus we 
have shown ([422|l and ((4231) . 

Furthermore, since e^^"^^!) = V, we have by the argument above (see (|4.30p and the argument 
after it) that (|4.24p holds. Moreover, Wtp S V, e~^^''^^ ^"^ip is analytic due to the formula e~^'^**'/^(/? = 
r(z)e^^'^°/^(/3 and analyticity of the factors on the r.h.s. . 

To prove ()4.25p we note that, due to the Dyson expansion, it is true for purely imaginary z and 
by analyticity of the l.h.s. for all z. 

Now, we prove ()4.26p . Denote by Yk{z)ip the truncated series on the r.h.s. of (I4.29p . with m < k. 
Let A G Tr{A)" ■ Choose a sequence of linear combinations of {Hp (8) Tip-valued) field operators 
An G 7r(y^) converging weakly to A. We use the defining property of the operator S = A~^/^ J, and 
the results proven above, to obtain the following relation for all (p E V: 

{JAT{z)no,^) = lim {JATkiz)na,^) 

k — >oo 

= lim lim {J AnTk{z)flo, ip) 

fe^oo ri— >oo 

= lim lim (rfc(z)*(A„)*f^o,Ai/V) (4-31) 
= lim (A*r!o,rfe(z)Ai/V, 
A*no,T{z)A'/^ip] 



We have used that (^„)*f7o E V C To, Ai/^i/. E V C J^o, and in ((OT|) that JA^/'^AnTkiz)no = 
Tk{z)* {An)*i^o- The latter relation follows from the facts that Tkiz) is affiliated with the von 
Neumann algebra Tr{A)" , and that {An)* leaves J-q invariant. It can also be verified directly, using 
the explicit actions of J and A^/^. This shows (|4.26p . 

To prove the last statement we introduce a new family of operators r^{z) :— e''Lo/2^-zL'- > /2 
related to the adjoint of T{z). First we prove that V C Dom(r'^(z)). To this end we note that, 
exactly as above, the formal expansion of r'^(z), which we denote by r'^(z), converges on elements 
of V and is entire as a function of z. Next, we observe that, since e^^°^'^'D = T>, the equation 
gzLo/2^ ^ Dom(e^^^'^*/^), proven above, implies that V C Dom(e^^-^*'^'/2). Let Vi := {i/' e H| = 
f{Lf))'ip for some / E C^(M)}. Then Vi is a core for e^^^"/^. Now we claim that for any tp E V 
and any ip E Vi, 

(e^^"/2^,e-^^"'/V) = (^,f^(^)^). (4.32) 

Indeed, the latter relation is true for z purely imaginary and it remains to be true for complex z 
by analyticity of both sides. The last relation and the fact that Vi is a core of the operator e^^"/^ 
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show that e-"^-^'"/2<^ G Dom(p/-^"/2) ^nd l^^e'^^l'^e-''^'''^ ^"^ip^ = {^,V^{z)^ . In other words, 

(/? € Dom(r^(z)) and r'^(z).^ = V^{z)Lp. This proves that V C Dom(r'^(z)) and that V^{z) is an 

entire operator-function. 

To prove it suffices to show that V{z)V D V. To prove the latter, let ■0 e I? C Dom(r'^(z)). 

Then e^^o/^e-^-^^'/s^ g Dom(e"^^'"/2e""^^''/2), and the image of this vector under e^-^*"/2g-^Lo/2 

is just ip. Hence tp S r(z)I' as required. 

This completes the proof of Lemma 14.21 ■ 
Remark: A sharper estimate can be obtained by using a standard argument to estimate the 

norms of the integrands (to keep notation simple we set ip = ^o)i 

II \hiOo{c^o"Hv)---ao"-{v)al^-{v)---al^^iv)y/^ (4.33) 

by using Wick's theorem, the expression for the imaginary-time two-point functions, in the same 
way as it is done, for instance, in [5], Thm IV. 3. 



5 Generator K and interpolating family A'(s) 

In this section we find an explicit form and some properties of the generator K of the one-parameter 
group U(t) introduced in the preceding section (cf. [21^), and of the family ^i'(s) which interpolates 
if to a selfadjoint operator. 

Let ojQ be the state of the algebra A fixed at the beginning of the Section |4] and let J and A be 
the Tomita-Takesaki modular conjugation and modular operator associated with the couple {A,u!o). 
We have the following standard relations: 

j7r(A) J 7r'(A), (5.1) 

JA^^'^T:{A)no = n{A*)no, A^/^I^q and A-**7r(A)A'* = 7r(K*(A)), where k is the automorphism 

of the algebra A defined in (|4.18p . The last three equations imply 

j7r(A)f7o -=7r(K*/2(A*))17o. (5.2) 

Finally, we recall that /? = max(/3i,/32) and that L^^^ is the self-adjoint operator defined as 
:= Lo + gV where V tt{v). 

Theorem 5.1 Assume that (A) and hold. The semigroup U{t), defined in is 

dijferentiable on the domain Dom(L^^)) fl TT(A)fl, and the generator K — —i-Sr U(t) is given on 
this domain by the expression 

K = Lo + giV-VU/^), (5.3) 

where V = Tr{v) and V/ = 7r'(7*(u)) with 7^* := ctg o (due to condition (A) the operator y'_ii2 
is well defined). Furthermore, the domain Dom(_L'-^^) n T:{A)fl is dense in Ti. 

Remark. The imaginary part of the generator K is not semi-bounded. Therefore, the group 
U{t), densely defined on 7r{A)^l, does not extend to a group of bounded operators. 

Proof of Theorem\5l\ Due to Condition (A), we have e-'^^o/^^.Wf^j, ^ ge-'^^'>/'^TT{v)no ^ 
g7r(CT^'^/^(w))f^o e 'D- It thus follows from Lemma[12]that L^'^^f^o E Dom(e"'3^'"/2). The latter fact 
implies that e'*'^*''ri is differentiable at i = and therefore fl g Dom(L'^^)). Now, let B £ 'n'iA) be 
such that BQ e Dom(L(^''). Taking into account equation (|4.3p we see that iJ^{B) and therefore, due 
to g3), also U{t)Bn, are differentiable in t at i = 0. Indeed, let i(e**-^'"Se-**-^''' - B)n ^F + G, 
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where F := e**^''' Bi(e-^*^*" - 1)^2 and G := i(e**^'" - l)Bn. Clearly, F ^ -iBL^^^fl and 
G iL^^^B^l. Now, differentiating equation (|4.2p we find 

KBn^ L'^^^Bn- BL'^'^^n. (5.4) 

Now we compute the last term on the r.h.s. of this expression. To this end we use the following 
relations: 

7r(A)f]o = Tr'{K'/^{A*))no (5.5) 

and, for z = it, 

e^^*''7r'(A)e-^^'" = 7:'{af{A)). (5.6) 

Eqn (|5.5p follows from relations (|5.ip and (|5.2p . Equality (|5.6p is proven by using the Kato- Trotter 
product formula. 

Now, we claim that Vip g T> 

{^,Be^^''\{v)n^) - ([7r'«oK^/2(^;))]*V,i3e^^'''^!o). (5.7) 

Observe that the vectors e^^* V = r(— 2z)e^^°iy9, G I?, are entire in z, by Lemma |4?2] Hence both 
sides of Eqn (|5.7p are entire in z. Therefore it suffices to prove (I5.7P for z = ii. Eqn (|5.7p with z = it 
follows from the relations (|5.5p and (|5.6p proven above. Thus (|5.7p is demonstrated. 

Now, let 1]' e-'^^'^'/^flo- Recall that L^^)^^' = ge-^'^*"/27r(i;)f7o. Then dEZ]) with z /3/2 
and the definition of the transformation 7^* imply that for all ip Cz V 

{iP,BL^^'>n') = g{^,Be-^^'-'^'^^{v)^la) 

= g{W{f'\vWi^,B^'). (5.8) 

This relation and the fact that D is a core for7r'(7*/2(-(;))* (which is a linear combination of creation 
and annihilation operators) show that BQ,' G Dom(7r'(7*/^(u)) and 

BL^'^^n = TT'{-f'''^{v))Bn. (5.9) 

Since the r.h.s. of (|5.9p is V^^^g^^i this equation together with (|5.4p implies (|5.3p . Finally, 

Dom(L(^^) n '!t{A)VL contains the set V and is therefore dense in Ti.. ■ 
Observe that 

K*|/3,=/3p=/3 = a^* and 7*|ft=/3p=/3 = id. (5.10) 

This implies that 

K\0^=p,=p = L, (5.11) 

where L := Lq + .9^ — 5^', with := 7r'(u), is the standard self-adjoint Liouville operator. In what 
follows we write X = Lq + gl , where 

i = v-vu,^. 

The operator K is non-self-adjoint for 5(3 7^ 0, and the perturbation / is not relatively bounded 
w.r.t. the unperturbed operator Lq. To study the evolution generated by K we use the family of 
operators 

:-Lo + 5(V^-K') , (5.12) 
where, recall, V — 7r(z;), and V!. :— ■k'{-^^[v)). This family interpolates between the operator K, 

K = if(_,/2) , (5.13) 
(see Eqn ()5.3p ) and the self-adjoint operators Ki^g-^ with real s. Under condition (A) on u, V^'j,-) and 
i^Tj-^-) are well defined on the dense domain Dom(L(^^) n ■k{A)^ for all s in the strip 

Se -.^ {t e <C\\luit\ < \ + e] , (5.14) 
for e > 0, and are strongly analytic there (recall that tt' is anti-linear). 
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6 Spectral Deformation of K and K(^s) 

Since the operator K is not self-adjoint it is not a simple matter to derive long-time properties of the 
dynamics e*^* from spectral properties oi K. As a result we bypass establishing the connection of 
the dynamics to the spectrum of K and instead connect it to certain spectral properties of a complex 
deformation, Kg, of this operator. To do this we use the interpolating family K(^s)e, which is the 
complex deformation of the family K(s) i (|5.12p . In this section we define complex deformations Kg 
and K(^gjg and in the next section we establish their spectral characteristics which are relevant for 
us. 

In order to carry out the spectral analysis of the operator K, which we begin in this section, we 
use the specifics of the Araki- Woods representation in an essential way. They were not used in an 
essential way for the developments up to this section. 

As a complex deformation we choose a combination of the complex dilation used in [81 and 
complex translation due to [18] (see [8], Section V.2 for a sketch of the relevant ideas). This complex 
deformation was used in [25j in the spectral analysis for a general class of Liouville type operators. 

First we define the group of dilations. Let Ud,5 be the second quantization of the one-parameter 
group 

Ud,S ■■ f{k) ^ e''/^f{e'k) 

of dilations on L^(R"). This group acts on creation and annihilation operators af{f) on the Fock 
space, 7i,., according to the rule 

Ud^5a*{f)U^} = a*{ud,sf), Ud^s^rj - f!,,. (6.1) 

We lift this group to the positive-temperature Hilbert space, (|2.33[) , according to the formula 

® Ud,s <E) Ud,-s ® Udj <E) Ud,^5- (6.2) 

Note that we could dilate each reservoir by a different amount. However, this does not give us any 
advantage, so to keep notation simple we use the same dilation parameter for both reservoirs. 

We record for future reference how the group Ud.s acts on the Liouville operator Lq and the 
positive-temperature photon number operator N :— 'Y^j=i -^j? where 

TVj := j [a^*j(/c)a£j(fc) +a,*j(fc)a^,j (fc)] Sk, (6.3) 

and the operators a^^ ^ (fc) were introduced after (|2.34p . We have (below we do not display the 
identity operators): 

Ud^sLrjUj^j = cosh{6)Lrj + sinh((5)Aj, (6.4) 
where A^- is the positive operator on the jth reservoir Hilbert space given by 

A,=/.(.)(a^,,(.)«.,(fc) + <,(.)a.„(.)).3,^ (,.3) 

and 

UdjN,U-j = N,. (6.6) 

Now we define a one-parameter group of translations. It can be defined as one-parameter group 
arising from transformations of the underlying physical space similarly to the dilation group. This 
is done in Appendix |B] We define here the translation group by means of the selfadjoint generator 
T Tj, where 

^ / [»lj(^)7a<'j(fc) + «*j(^)7arj(^)] d^k. (6.7) 
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The operator 7 = i{k • V + V • A;), with k — k/\k\, is a symmetric, but not a self-adjoint operator. 
Nevertheless, the operators Tj,j = 1,2, are self-adjoint [25j. Thus the operator T is self-adjoint as 
well. We define the one-parameter group of translations as 

Ut.r ■■= lp®lp® e'^^. (6.8) 

Eqns. (|6.7p - (|6.8p imply the following expressions for the action of this group on the Liouville 
operators: 

Ut,rLrjU[::^ = Lrj + tNj, (6.9) 

and Ut^rNjU^^ = Nj. Observe that neither the dilation nor the translation group affects the particle 
vectors. 

Now we want to apply the product of these transformations to the full operator K = Lq + gl. 
Since the dilation and translation transformations do not commute we have to choose the order 
in which we apply them. Since the operator A = Aj is not analytic under the translations 
while the operator N is analytic under dilations we apply first the translation and then the dilation 
transformation. We define the combined translation-dilation transformation as 



Ue = UM.r (6.10) 

where 9 — (^, r). Note that Ug leaves the finite-particle space J-q, as well as Dom(A) n Dom(A^) 
invariant, for 9 €M? . 

In what follows we will use the notation 16*1 = (|(5|, |t|), Im0 = (Im(5, Imr), and similarly for Re6', 
and 

Im^ > Im(5 > A Imr > 0. (6.11) 

Now we are ready to define a complex deformation of the operator K . On the set Dom(A) n 
Dom(iV) we define for 9 

Ke := UgKUgK (6.12) 

Recalling the decomposition K — Lo+gl, where Lq Lp + Lr, Lr '■— ^rj and / = V — V^^^^^ 

we have 

Ke^Lo,e+9le, (6.13) 
where the families Lo g and Ig are defined accordingly. Due to Eqns. (|6.4|) . (|6.6p and (|6.9|) we have: 

Lo g = Lp + cosh{S)Lr + sinh((5)A -I- tN, (6.14) 

where 9 — (5, t), and A = j^i ^j- explicit expression for the family Ig is given in Appendix 
|R2](see Eqns (IbTs]) and ^^J^i). 

Similarly, we define the family K(^s)0 UgK(^s)Ug^ (recall (|5.12p ). 

The operator families above are well defined for real 9. Our task is to define them as analytic 
families on the strips 

= {61 e C^IO < ±^116* < 6*0} (6.15) 

where 9o — {So, To) > is the same as in Condition (B). Recall that the inequality ±.Im9 < 9o is 
equivalent to the following inequalities: ±Im5 < Sq and ztlmr < tq. (The fact that analyticity in a 
neighbourhood of a fixed 9 g S*^ implies analyticity in the corresponding strip in which Re0 is not 
constraint follows from the explicit formulas (|6.14p . (|B.2.5P and (|B.2.7p .) The analytic continuations 
(if they exist) are denoted by the same symbols. 

We define the family Kg for 9 e {9 & C2||Im6'| < 9o} by the exphcit expressions (pTT^ . ((Oi)) . 
(|B.2.5P and (|B.2.7p . Clearly, Dom(A) n Dom(iV) C Dom(Lo6i) and on this domain the family Log 
is manifestly strongly analytic in 9 £ {9 £ C^||Im6'| < 9o}- It is shown in Appendix [B] that for 
\lm9\ < 9o we have Dom(A^/^) C Dom(/e) and Igf is analytic V/ S Dom(A^/2). Here Condition 
(B) of Section [3] is used. Hence the family Kg for 9 G {9 & C^||lm0| < ^o} is bounded from 
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Dom(A)nDom(7V) to V. (and Kef is analytic in 6* £ {6* £ C^||Ini6l| < Oo}, V/ e Dom(A)nDom(iV)). 
Moreover, for jlm 0\ > the operators Ke are closable on the domain Dom(A) n Doni(iV), since 
their adjoints are defined on this dense domain. We denote their closure by the same symbols. 

However, {Kg \ IvciO < Oq} is not an analytic family in the sense of Kato. The problem here is 
the lack of coercivity - the perturbation / is not bounded relatively to the unperturbed operator 
Lq. To compensate for this we have chosen the deformation Ue in such a way that the operator 
Mg := ImLo,e is coercive for lm9 > , i.e., the perturbation Ig as well as ReLp^e are bounded 
relative to this operator. The problem here is that Mg — > as Im6' — > so we have to proceed 
carefully. 

Everything said about Kg applies also to the family K(^s)e- 

The next result gives some analyticity properties and some global spectral properties of Kf^g-^g. 

Theorem 6.1 ([25]) Assume that Condition (B) holds and let 6^ — {5q^t^) he as in that condition. 
Take an 



a > 



9 



2 



sin(Im(5) 



where Cq := C(l + + P2 ^^^); '^^^ where C is a constant depending only on tanSo. Then we 

have: 

(i) {z G C|Imz < —a} C p[Ki^g)g) (resolvent set) if either s G M and 9 £ , or if s £ Se and 

(ii) Let Ca.b be the truncated cone 

Ca,b := {z e C I Imz > -f, |Rez| < 2[{smb)-'^ + a/A]{lm z + a) + \\Lp\\ + l} . 

Take s £ Se, £ S*^, and take a as in S6.16\) . Then a^Kg) C Ca.imS, and for z £ C\CQ,iin5 
we have 

\\{Kg-z)-^\\ < (dist(z,a,Im5))"'. (6.17) 

(Hi) The family K(^g-^g, s £ Se, £ S*^, is analytic of type A (in the sense of Kato) in 6 £ 5^, for 
all s £ Se, and in s £ S^, for all 9 £ 5^; 

(iv) Let s £ M. For any u and v which are Ug-analytic in a strip G C^|0 < ImB < for some 
Oi = ((5i,To), Si £ [0, min{7r/3, 0o})> the following relation holds: 

{u, (K^g) - z)-\) - {u^, iK^g)g - z)-\g) , (6.18) 

where ug = Ugu, etc., Imz < —a and < lm6 < 9i/2. 

Proof. Statements (i), (iii) and (iv) are special cases of Theorem 5.1 in [25j . with the exception of 
the assertion about analyticity of s i— > K(^g^g in (iii). This assertion is easily proven by noticing that 
ds{Ki^s-^g — z)~^ = —{K(^g-^g — z)^^ {dsVg){K (^g)g ~ z)^^ . Statement (ii) is the content of Proposition 
5.2 in ESI. ■ 



7 Spectral Analysis of Kq 

In what follows, given a self-adjoint operator A and a g R we use the notation XA<a for the spectral 
projection of A associated to the set {A G K |A < a} and similarly for XA>a and XA=a, etc. Fix a 9 
satisfying < Im6' < ^o- 
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Theorem 7.1 Assume Conditions (A) - (D). Let 9q be the same as in Condition (B), let gi be as 
m Theorem\3Ji and let Bq > Ime = {S' , t') > 0. If < \g\ < gi and t' > \g\'^+" then 

(a) is an isolated and simple eigenvalue of Kg; 

(b) (7{Kg)\{Q} C {z e C+|Imz > min(c52, \t')}, 
for some c > 0, independent of 0. 

Theorem lT.ll is proven at the end of this section. Together with Theorem l6.1l it gives the following 
picture for the spectrum of Kg. 




Spectrum of Kg for Im6' > 



The motion of resonances bifurcating out of the eigenvalues of Lq is governed, to second order in 
the coupling constant by level shift operators, see |25l 124] for a discussion closest to the situation 
at hand. Let e be an eigenvalue of Lp and let Ag be the level shift operator acting on Y^a,r^XLp=e^ 
defined by 

Ae:--Pe/(io-e + iO)-i/Pe, (7.1) 



where Pe = Xip=e ® Xl^=q- The notation +iO in (|7.ip stands for the limit of as e J, 0. The 
following result summarizes properties of the level shift operators which are essential for the proof 
of Theorem O 

Theorem 7.2 Assume Condition (C). The level shift operators Ag, |7.-?[), satisfy 



(T(Ae) c C+ (7.2) 

.(A.)nM={J^^ lltl (7.3) 

dimKer(Ae=o) = 1 • (7.4) 

Furthermore, there is a "fo > which does not depend on the inverse temperatures, s.t. ImAg > 70 ; 
for all e 0. Moreover, if Condition (D) is satisfied, then there is a 5q > Q which does not depend 
on the inverse temperatures, s.t. Ini((T(Ao)\{0}) > 5{). 

We prove this theorem is in Section [TUl 

Proof of Theorem \7.1\ That is an eigenvalue of Kg follows readily from the equations Kil = 
and the fact that Q is J/g-analytic in the strip (|6.15p . as we show in Lemma [731 below. So we have 

Kgilg = 0, (7.5) 

where fig := Ugil.. 

Lemma 7.3 il is Ug-analytic, for 9 — {S,t) G , see i6.15]) . 



21 



Proof of Lemma \77^ This follows from the Dyson series expansion ()4.29p given in Lemma [4.21 
and the analyticity condition (B). ■ 
Let po G (0, (7/2), where a is given in (j3.9p . and consider the half-space 

5 = {z e C|lmz < iposin(Im(5)/9o} ■ (7.6) 

We decompose this region into the strips Se ~ {z E S\ |Rez — e| < po}: where e e a{Lp), and the 
complement. The following result is a special case of Theorem 6.1 of [25]. 

Theorem 7.4 Assume that condition (B) holds and that there are constants 70 > and 5o > Q 
satisfying ImA^ > 70 for all Q ^ e <E <j{Lp) and Im((T(Ao)\{0}) > 5q. Take < I5I < y^go and let 
a — {11 — l/2)/(/i + 1/2), where /i > 1/2 is given in Condition (B). Then 

1. Wehavea{Ke)nS C Uee^(L^) ^e- 

2. Choose po = There is a C > s.t. if < \g\ < C{jo)^^°', then, for all e^O, 

(j{Kg) n 5e C {z e C I Imz > ig2^o}, (7.7) 

3. Choose po = IffP"^"- There is a C > s.t. ifO< \g\ < C min[((5o)i/", (r')^] then 

a{Kg) n 5o C {zo} U {z G C I Imz > i mm{g'^do, t')}, (7.8) 
where zq is a simple isolated eigenvalue of Kg, satisfying |zo| = 0(|(7p^"). 
For a coupling constant satisfying 



< l^l < min 



(5o)'/",(7o)'/",(^o)'/",(r')^l (7.9) 



all three parts of Theorem 17.41 apply. Thus Kg has a simple isolated eigenvalue {zq} in a neigh- 
bourhood 0(|gp+") of the origin, and a{Kg)\{zo} C {z G C+ | Imz > mm{cg^, ^t')}, where 
c — i min(7o, (5o). In order to complete the proof of Theorem 17.11 we only need to remark that 
Zq = since zero is an eigenvalue of Kg, see (17. 5p . (One can also give a dynamical argument to 
prove that zq = 0, see the remark after (I8.12[) .) ■ 



8 Resolvent representation and pole approximation 

In order to study the long-time behaviour of the evolution U (t) — e'^* (defined on the domain 
7r{A)^l by (|4.2p ') we relate it to an object which we understand relatively well, namely the resolvent 
{Kg — z)"^ of the deformation Kg, defined in Section [HI The main result of this section is 

Theorem 8.1 Assume that Conditions (A), (B) and J^. J[ j hold. Let (p and he Ug-analytic vectors, 
and let — tt{A)Q for some A ^ A. We have the following representation 

{<P, e^^**) ^^l-j{cjy^, {Kg ~ zY^^g) e^^'dz , (8.1) 
r 

where = Ug'i' (similarly for (p or any other vector), 9 G S*^, and F is the path 

F := {z = A - ir'/a, | A| < C} U {z = A - i2T' /?, + iAr'/SC, A > C} 

^{z = \-i2T' /i-iXr' /iC,X<-C} (8.2) 

for a sufficiently large constant C . The integral on the r.h.s. of hS.l]) is well defined in virtue of 
Theorem \6.1[ ^6.17^ , and the estimate |e*^*| < e"'*''^ /'^'^ on the infinite branches ofT. 
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Proof. In the proof below, the vectors (j) and 4" are as in the theorem. To prove the equahty in 
()8.ip we use the family K(s) of operators defined in ()5.12p - ()5.13p (here we use Conditions (A) and 
(gH)). Note that Q e Dom(e''^) for any ly > 0, as follows from the relation G Dom(e''^e-'^^'''/2) 
which is shown in the same way as the relation Qq G Dom(e~''^**'/^), see Lemma 

Next, we define the operator e'^t")* as follows: e*^('')*ri is given by a Dyson expansion, where the 
part —gVg of K/^g-j is treated as a perturbation. The fact that the Dyson series converges is easily 
seen from the relation il e Dom(e'^^), > 0, shown above. Moreover, it is clear that this series 
defines a vector which is analytic in s € S'e, i.e., s i-^ e*^<=''ri is analytic for s e S,.. We define the 
action of e'^(=)* on vectors Tr{A)n, A G A (which form a dense set), by 

e'^M*Tr{A)n = cr*(7r(A))e*^<=)*f7. (8.3) 

Consequently, the map s i-^ e*^<=)*7r(^)f7 is analytic for s e S^. For s = ^ this definition gives 

Since -f^'(s) is self-adjoint for s real we derive from Stone's formula 

(0, e^^c^)**) ^ __L ^ (0, (K^,^ - z)-H) e"-'dz, s G M. (8.4) 

K-i 

Next, using e*^* = ^ a(Rez) ^'^* ^^"^ integrating by parts we can represent the r.h.s. of (|8.4p as 



Now we perform the spectral deformation. Theorem 16.11 (iv), to obtain for 9i/2 > ImO > (where 
6*1 is given in Theorem 16. 11 (iv); here we use Conditions (B)) 

RHS^^-^^^^ J {c^,{Kf^,)e~z)-He)e''dz. (8.5) 

The integral converges since due to ()6.17p we have ||(i4r(s)e — ^:)~"|| < C„(Rez)^" for z e R — ie, 
where {x) := VT + a?. 

In (|8.5p we deform the contour of integration from R — z to F which is fine due to (|6.17p , and we 
integrate by parts in the opposite direction of above, to obtain 

(0, e'^(=)*M') = j (<%, (K(,)e - zyHe) e'^'dz . (8.6) 

r 

Both sides of this expression are well defined and analytic for s £ (see Theorem 16.11 ()6.17p . and 
after equation (|8.3p ). Since they are equal for real s they are equal for all s in their domain of 
analyticity. Taking s = — | in this formula gives (|8.ip . | 
It is shown in Section [7] that the operator family Kg has a simple isolated eigenvalue at and 
the rest of its spectrum is located in a truncated cone in {z G C~''|Imz > ^t'}, where r' = Imr. In 
the integral on the r.h.s. of formula (|8.ip we deform the contour of integration to 

r' ■.= {z^X + iir', |A| < C} U {z = A + iAr'/SC, A > C} U {z = A - iXr' /3C\ A < -C} (8.7) 

where C is sufficiently large. Picking up the residue from the simple eigenvalue of Kg we derive 
from dHl]) 

(0, e'^**) = (<%, Pg^g) + 0(110^11 II II 6--'*/=^) , (8.8) 
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where Pe = ^ §{^9 — z)~^dz is the eigenprojection of Kq corresponding to the simple and isolated 
eigenvalue and the remainder bound is coming from the term ^ /p/ (</%, [Kg — z)~^'^g) e^^^dz. 
The contour integral is over a small circle around the origin and the path V is defined in (|8.7p . 
This is the only place where we use that is an isolated eigenvalue of Kg ~ the fact we show using 
complex translation in addition to complex dilation. 

Pg is a rank one projection which is analytic in 6 G 5^. One proves using a standard argument 
that it satisfies KgPg = PeKg = 0. Hence, Pg can be written as 

Pg = \ng) (r!l| , (8.9) 

where fig and fJi are zero eigenvectors of Kg and its adjoint operator 

{Kg)* ^{K*)g=:K^, (8.10) 

i.e. we have 

KgQg = and K^fl^ = 0, (8.11) 

with the normalization ilg) — 1. Since Pg and fig are analytic in 9 G S*^, then so is fl^ in the 
variable G S'g^. (The possibility of the normalization, {fl^,flg) ^ 0, follows also from results of 
Section [11] Analyticity of fl'^ in G S"^ can also be shown directly by using the analyticity and 
spectral properties of K^, see Section [7|)- Equation (|8.9p implies 

{<j>g,Pg^g) = n){n*^,^g) . (8.12) 

The resonance vector fl* appearing in (|2.32p and (|4.13p is defined as 

{^,n*) = {<i>_^,n*g) . (8.13) 

Remark. We present here another proof of the relation zq ~ 0, where zq is the simple isolated 
eigenvalue of Kg given in Theorem 17.41 Starting with the information on the spectrum of Kg given 
in Theorem 17.41 and proceeding with a contour deformation as above we find that (0, e*-f^*V) = 
e^^^n* {(j>g, Pg^pg) + 0(e(|i™ ^"l-^'/^)*) instead of ([g^gjl . Applying this formula to (j) = = fl and using 
g3) we obtain (O, ^2) = e*^«* (17, rj) +0(e(|i™^«l-^'/3)t). Since |zo| = 0(|5p+") « t' the error term 
tends to zero as t — *■ oo, so by taking limr^oo Jq dt on both sides we see that zq must be zero. 



9 Proof of Theorem [33] 

Let Tp be any wo-normal state on A. For the following reasoning, we may assume without loss of 
generality that ip{A) = (fJ^, 7r(A)ri^), for some fl^f, G H. Since fl is cyclic for Tr{Ay , ip can be 
approximated as follows. For any e > there is a Q G i^iA)' s.t., for all A £ A, 

\^{A)~ {Qfl,n{A)Qfl)\<e\\A\\. (9.1) 

Applying this to ijj^{A) — {fl^^a^{Ti{A))Vljp) , pulling Q through (T*(7r(j4)) and taking into account 
(^301) . U{t) = e"-^ and U{i)fl = 0, we obtain 

\i^\A) - (Q*Q17,e**^^(A)f2) I < eUl (9.2) 

uniformly in i G K. In order to examine the long time behaviour of {Q*Qfl, e'*'^7r(yl)r2), via (|8.8p . we 
first approximate the vector Q*Qfl G W by a family of C/g-analytic vectors s.t. ||xe — Q*Qfl\\ < e. 
We have 

\{Q*Qn-Xe.e^'''7TiA)n)\<e\\e^'''TriAM^e\\a\niA)M<e\\A\\. (9.3) 
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It follows from ((921), (EH) and Kl2\i that 

\^\A)-{x.,n){n^,{7T{A)n)o)\<2e\\A\\+C\\x,^^\\ \\{7r{A)n)e\\e--''^'. (9.4) 

Since {Q*Qil, fi) ^ 1 as e ~> 0, we have {xe, ri) = 1 + o(e"), where o{e'~') denotes a quantity that 
vanishes in the limit e ^ 0. Thus 

|V*(A) - (x., 1^) {nt {7riA)n)e) \ > \^P\A) {Qt (^(A)f]),) | - o(6°) | {^t (7r(A)17),) |. 

Combining this estimate with (|9.4p we arrive at 

|V'*(A)-(f71 (7r(A)17),)| 

< 0(6°) (||A|| + I {n^, i7:iA)n)e) |) + C||x,,^|| ||(7^(A)^7),||e--'*/^ (9.5) 

where o{e^) is independent of A and 6*. Let rj be the state on A given by 

r;(A) (^(A)f}),) . (9.6) 

Note that ri{A) is independent of the deformation parameter 9, if < Im^ < Oq, and < g < go 
{go depends on sin(Imi5), see the equation after (|3.9p ). Taking in (|9.5p first t ^ oo and then e ^ 
yields 

lim V*(^) =?7(^), (9.7) 

t — *oo 

Equation (|9.7|) shows that | (^fi^, {Tr{A)Q)g'^ \ < \\A\\. We re-inject this inequality into r.h.s. of (|9.5p 
to arrive at 

\i;\A)-n(A)\ , , 

lim sup ''^ ^ \ ^' =0, 9.8 

*-°°A6^ infill 

where 

^= {A e I 7r(A)f7 is L/e-analytic for \B\ < 9o}, (9.9) 
and where 1 1 1 • 1 1 1 is the norm on A defined by 

--\\A\\+ sup \\{7T{A)n)g\\, for AG A (9.10) 
\o\<eo 



Observe that 1 G A, and that the normalization (^il^^flgj = 1 implies 77(1) = 1. is a linear 

subspace of A, but not an algebra. 

We show in Appendix [Cl Proposition IC.li that Ai (defined in (|2.17p ) is strongly dense in A 
(defined in (|2.7p ). and that any A G Ai has the property that Tr{A)n is [/e-analytic, for 61 in a 
neighbourhood of 6* = 0. Hence Ai A A, and therefore A is strongly dense in A. 

We have thus shown that for any wp-normal state -0, i/;* — > 77 as t — > 00, where the convergence 
is understood in the ||| • |||-topology of linear functionals on A. 

It is clear from (19. 8p that |77(A)| < \\A\\, for A £ A, hence extends to a bounded positive linear 
functional on the Banach space of observables 

Ao ■■= II • ||-closure of A, (9.11) 

normalized as 77(1) = 1. Standard perturbation theory shows that 77(A), A £ A, is analytic in the 
coupling constant g. 

Observe that we can rewrite the state 77(A) also in the form 

T^iA) = Tr {TriA)gPe) (9.12) 

where TT{A)g := Ug'K{A)Ug^ . Formally one can undo the rotation in (|9.6p to obtain (|4.14p with 
n* := Uj^Qj. However, in the non-equilibrium situation fl* ^ W. The set Ai is exactly the set on 
which (|4.14p makes sense. Thus we gave a rigorous meaning to (|4.14p and the NESS vector Q* . 
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10 Proof of Theorem [72] 



Our task is to show that the spectrum of (|7.ip lies in the upper complex half plane {Im z > 0} if e 7^ 0; 
and that it has a simple eigenvalue at zero and all the other eigenvalues lie in the upper complex 
half plane if e = 0. While this analysis is standard in the case when / is a selfadjoint operator (then 
the imaginary part of (|7.ip is just PeIS{Lo — e)IPe, manifestly a non-negative operator; see e.g. [8]), 
it needs some more thought in our case, where the interaction is non-selfadjoint. Let 

V,=niv,), and V^^n'{v,). 

The main ingredient of the proof is 

Proposition 10.1 Assume Condition (A). We have 



Hp) A,e (e-f"'"-'^ ® tp 



(10.1) 



where, setting TZ := {Lq — e + iO) 



A., 



^p,{v,-v;)n{v,-v;)p. 



(10.2) 



Notice that UP. 1\) shows that the spectrum of Ag is independent of Pp . 



The importance of (|10.ip is that it relates Ae to the operators Aje whose spectral characteristics 
are well known. Indeed, Aje are the level shift operators corresponding to the reservoir j coupled to 
the particle system, studied in [51 [Mj ■ 

Before proceeding to the proof we examine some consequences of this proposition. We assume 
Conditions (C) and (D) in addition to Condition (A). 

The case e ^ 0. Let us assume that the nonzero eigenvalues of Lp are simple, i.e. Ei — Ej = 



E„ 



En ^ i ^ m, j = n. For a treatment of the more general case where Ei — Ej = E„ 



En 



for («,j) ^ {m,n), with simple Ej, we refer to |26| . Since Pe is of rank one Ag is just a complex 
number, namely the sum of Aie + A2e (the dependence on f3i, P2 disappears). Under condition (|3.5p . 
one has ImAg > 70, where 70 is a strictly positive constant which is independent of the inverse 
temperatures, see [5^). This shows (|7.2p and (|7.3p for e ^ 0. 

The case e = 0. Zero is necessarily a degenerate eigenvalue of Lq, so the above reasoning does 
not apply. In particular, Aio and A20 do not commute. It is shown in |24[ 126] that 



Ajo = ilmAjo iL 



jO, 



(10.3) 



where Im A_,o '■— ^{-^jo — ^jo)- We use here implicitly Condition (C) on the non-degeneracy of Hp; 
if the small system has degenerate energy levels then AjQ are not purely imaginary [24]. One shows 
as in [8l[26ll24] that Tjo are real matrices having strictly negative off-diagonal entries, (rjo)mji < 0, 
for m ^ n, and satisfying 

Tj^nf^^^O, (10.4) 
where ^l^p'^ is the particle Gibbs state at temperature [3j. Hence, since 



we see that 



i=i,2 



pj V tTe-f'.Hj2 r.of^p ' -0, 
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where Tq := —iAq. Thus, ilp^''^ is an eigenvector of the real matrix Tq with eigenvalue zero. Notice 
that the vector uf'''^ has strictly positive components, [tre~''^^^/^]~^e~'^''^"/^, in the orthonormal 
basis {ifin'^^n} of RanxLp=o (where Hpipn = En^Pn, \\Vn\\ = !)■ Moreover, the off-diagonal elements 
of the real matrix (which is not symmetric for Pi (^2) To are given by 

(roU„= J2 [(e(-^''+^^-)^-/^®]ip)r,o(e-(-^''+ft)^-/2 0]ip)]^^ 

j=l,2 

for m > n (and similarly for m < n, see also [SJIIS]). Hence Condition (D) implies that {To)m.n < 
0. A standard Perron- Frobenius argument shows that zero is a simple eigenvalue of Fq, and that 
cr(ro)\{0} C C+. This shows equations (|7.3p - (|7.4p for e = 0. It is shown in [8] that the gap at 
the bottom of the spectrum of TjQ, j = 1,2 has a lower bound which is independent of the inverse 
temperatures. 

We now prove the existence of 60, assuming that Condition (E) is satisfied. If dim Tip = 2 then 
one eigenvalue of Aq is zero and the other equals the trace of Aq. Expressions (|10.ip and (|10.3[) show 
that Tr(Ao) = «[Tr(rio) + Tr(r2o)]. Thus the spectral gap of Aq is the sum of the gaps of Foi and 
ro2, which have lower bounds uniform in the inverse temperatures. 

Next take dim Hp > 3. For 6(3 = \Pi — ^2] =0 the matrix Aq has the same spectrum as J2j ^jo, 
see (|10.ip . An application of the minimax principle demonstrates that the spectral gap of the latter 
operator has to be at least as large as the maximum of the gaps of Tjo, j = 1,2. For small values of 
SP (c.f. (|3.6p ). the existence of So follows by perturbation theory. 

Finally we consider the case where Sp and Pi, P2 are large (see p.6p ). Let us take P2 = Pi + Sp. 
As is easily seen from (jlO.ip we have 

a(Ao) = la ( [Fio + {e^P""'^ ® t)T2o{e~^P"'''^ ® 1)] ) . (10.6) 

Using the explicit expression for the matrix elements of in the basis Lpj ® ipj (which can be 
read off of relation (jlO.Sp for off-diagonal terms, and is easy to obtain for the diagonal ones), one 
verifies that the matrix (e'^'^^p^^^ lI)F2o(e~''''^p/^ (g) 1) converges to a lower-triangular matrix Q{Pi), 
in the limit (5/3 — > cxd (uniformly in Pi), and furthermore, that Q{Pi) D as Pi ^ 00, where D 
is a diagonal matrix with non-negative entries. The minimax principle implies that all but one 
eigenvalues of Fio + D are greater than, or equal to the gap of Fio. From perturbation theory we 
know that for &P and Pi sufficiently large (independently of each other), all but one eigenvalues of 
the operator Fio -I- [e^P^p/^ (g) ll)F2o(e~'''^^p/^ g) 1) must have real part greater than, or equal to half 
of the gap of Fiq. The existence of Sq now follows from (|10.6p . ■ 
Proof of Proposition [TUA\ Let V — ^ and V — J2^=i ■ definition of the 

operator I , I := V — VLi/2J and the relation V^^^2 — e^^^V'e^^^^ we have 

where L := SPpLp + SPiLri + 6P2Lr2 with SPp = Pp-P and SPj = Pj - p. Using that P^e-^l^"^"^^ = 
P^e^^P'''^^^ we decompose 

Ae = PeVnVPe + PeV'nV'Pe 

-PeU7^e■^/2y'p^g-5/3pe/2 „ p^v'e-^^'^nVPee^'^''''/'^. (10.7) 
Notice that V and TZ commute with Hp (g) g^l^pHpf^^ Using this and the relation 

(Ip (g) e-^/^pifp/a^p^ ^ g-5/3pe/2^g5/3pHp/2 ^ ^^^p^^ 
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we obtain 

Next, using (|10.8p again, we find 
P.VTZe'^'^V'Pee-^P-^'^ 

-S(i,H,/2 ^ p^y^^{Sp^L^^+S0,L^,)/2y,p^ (^^S0,H,/2 ^ 

Treating the other two terms in (|10.7p in a similar way, we arrive at 



-Ae = [e-^p"""/^ ® iipj p^ivnv + v'nv 

We now examine the term in [•••]. Write F = Vi + T^, l^' = V/ + Vj, where, recall, V, = 7r(uj), 
1/' = -K^Vj). Notice that we have {\®\®Pi-i®\2)y\{\®\®PQ.®'^2)=^, and similarly for 
V2, from which it follows that the expression Pe[- • • ]Pe in (|10.9|) splits into a sum 



Pe WjTZVj + V'nV^ - -^/.7^e('5/3li-l+'5&L.2)/2y^/ _ y^'j^^i-SPiL^i-Sp2L.2)/2y^ 



i=l,2 



P. 



We consider the j — 1 term. Using that Lr2 commutes with Vi, V{ and that Lr2Pe ~ 0, we see that 

Pe^l7^e'''^l^•■l+'''^2-^'■^y/Pe = Pe^47^e'^'^l^•■l/2y/Pe ^ PeVine^^^^°/'^e^^'^^^^''^v[Pe 

= (e^^i^^'/^ ® ]lp) P.ViTZe'^'^'^^-'^^^VlPe (e^'P^"-'^ ® Hp) . 

All other terms in p0.9p for j — 1, as well as the terms for j — 2, are treated similarly and one 
arrives at 



-Ae = (^e-'^^"^/'^ (g) Hp) ^ (^e^^Hp/^ <g, ]l) PelVjUVj + V^UV^ 

-VjTle'^^^'^^-'^^^V; - y;7^e-*ft(^"-'=)/VJ]Pe (e^ft-^-/^ ® Hp) (e/3p^p/2 ^ , 
where we used SPp - 6Pj ^ l3p- 13-j. Hence, -Ae = RHS ^JT^ + {e-^p"p ® ]lp)i?(e'^p^f (g) Ip), where 
P = ^ (e'^^^''/2® ]l)Pe[y,7^(l-e*^^(^°-^)/2)y/ 

+y;7^(l - e-*'^^(^°-^)/')Fj]Pe (e^ft^^/^ ip) . 
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Since Lq implements the free dynamics, we have that e^^^°Vje '^^^ commutes with Vj , for z €E C. 
Using this, writing ^Lo-e±io — ~ ~ Jo dse^^ " ^' and using that Pg-^o = ePe, we see that 

1 _ g5ft(io-e)/2 1 _ -SPj{Lo-e)/2 

PV V'P — — PT/' - VP 

Lo - e - iO ■' ■' Lo- e + iQ 

Consequently, R = 0. This concludes the proof of Proposition llO.il | 



11 Perturbation Theory for NESS 

In this section we develop a perturbation theory for the NESS r]{A) := (^^, '!T{A)gflg^ . The vectors 
n^- and fig are the zero eigenvectors of the operators and Kg respectively, see ()8.10p and ()8.1ip . 
We derive perturbation expansions for and fig, see (|ll.lip and (|11.15p below, using the Feshbach 
maps introduced in [6l [7] , and extended in [5] . We review the definitions and some properties of 
these maps referring the reader to [71 [5] for more detail. For simplicity we present here the original 
version, [6l[7], though the refined one, [5|, the smooth Feshbach map, is easier to use from a technical 
point of view. 

Let X be a Banach space and let P be a projection on X. Define P := I P and let H-p PHP 
and Rp{H) := PH^^P if H-p is invcrtible on RanP. We define the Feshbach map Fp by the relation 
Fp{H) -.^ P{H - HRp{H)H) P on the domain 

Dom(Pp) = 

{H -.X X\Hp is invcrtible, RanP C Dom(iJ ) , RanP-p (i? ) C Dom(Pi/P)}. 
A key property of the maps Fp is given in the following statement proven in [7]: 
Theorem 11.1 (Isospectrality Theorem) ("ij G a{H) G a{Fp{H)), 

(ii) Hij = ^ Fp{H)ip = with ip = Pi, ("^") and V' = (1 - Rp{H)H)ip ("<^"). 
Let Pep be defined as 

Pep XL^=e «) XMg<p, (H-l) 

where XLp=e is the eigenprojection for the operator Lp corresponding to an eigenvalue e G o'{Lp) 
and XMg<p is the spectral projection for the self-adjoint operator Mg := ImLo.e corresponding to 
the spectral interval [0,p] (remember that Mg is a positive operator). 
The following result is proven in [25], Lemma 6.3. 

Lemma 11.2 Assume Condition (B). Take po G (0, (t/2) and let \g\ < ^/pogo, where go is given 

after ()3.9|) . If z E Se then Kgz '■— Kg — z G Dom(Pp^ ), and the operator Kg]} := Fp^ (Kgz) 
acting on RanPgpg is of the form 

Kf] = (e - z)l + Lrg + g^Ae + 0{e{g, po)). (11.2) 

The remainder is estimated in operator norm, \\0{e{g, po))\\ < Ct{g^i^, po), where g^^^ = 
|g| maxj sup|g|<g^ IjGj 115^^.^1/2,9, with a C independent of 5[3j, 6, and where we have set 

e{9, P) \g\p' + \9\'p-"'- + \g? p"'-' ■ (11-3) 
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To unify the following analysis we write ilf for either fig or fii. Correspondingly, I"^ , Lf^gn. and 
Kf stand for either Ig or /i, for either L^q or Lp^, and for either Kg or K^, respectively. We use 
the shorthand Pq = ^Opo i?o(^) •= Pa^-p^Poi where Ap := PAP. We will assume that 



1/2' 

Theorem 1 1 1 . 1 1 and Lemma [11.21 imply that Kf £ Dom(Pp„), that 

Fp„{K*)Pon* ^0, 
and that the original eigenvector il.f can be reconstructed as 



nf = 



l-gRo{K*)I* 



We expand Ro{Kf) in this expression into a Neumann series, 



(11.4) 



(11.5) 



(11.6) 



(11.7) 



~ 1 / 2 

for any iV > 1, provided that 0{{g^i^pQ )) = Og(l). The remainder term in (|11.7p is obtained by 
using a standard estimate on the Ai"th term of the convergent Neumann series. Indeed, writing 



N 



Po{Poe*)lf 
and using the estimates 

|5| (Me+po)-'/'/f(Me + Po)-'/' 



(11. 



=-Po(Mg + po) ' Ig[Mg + poj ' 

^ne*Po 



1 N 



[Mg + Po) 



1/2 



< 



1/2 



(see also Lemma 5.3 of [25j) and 



Mo+po - 



< C, we obtain 



N 



Po^ 



# 



^ -\ 2, ~1 2^N 1/ 

<P0 idsfiPo ) Po 



(11.9) 



(11.10) 



(11.11) 



Observe that since g^ffPo^ = Og(l) we have 

- E 9''i-MLo0*)lfrPont + oig). 

n=0,l 

Let := Fp„{K*). As in (flT^ it can be written as 

where K*^^^ L^e* + 5^A# with A# either Aq or Aq, where Aq is given by (fTTj) with e = 0. We 
assume that Sf3 varies in the set \Sp\ < c, for some c > 0, so that g,^ can be replaced by g. Take 



Po = 9 



with a = then Lemma 111 .21 gives 



W 



0(5'+"). 



(11.12) 
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By Theorem 1 11.11 

K*^^^Ponf^O. (11.13) 
Let Qf :— Xx#(i>-o ~ Xa#=o ® XLr=Qi where Xa#=o is the Riesz projection onto the kernel of A'^. 

The operator Kq is normal with the simple eigenvalue separated from the rest of the spectrum 
by a gap > cmin(r', g^) for some c > 0. The operator A' ^'^^^ has also the simple eigenvalue 0, which, 
by the Kato-Rellich theorem, (I11.12[) and the assumption (|11.4[) . is separated from the rest of the 
spectrum by a gap > c'min(T', g^) for some c' > 0. Hence, we conclude that X'^'^^ G Dom(i^„#). 
Therefore, by Theorem 1 11.1) 

Ponf = (i - R-^{K*(^^))w*) Qfnf. (11.14) 

Since \\R-^{K'^^^'>)\\ < C[min(T', y^)]^^ we have the absolutely convergent perturbation expansion 

Qo 

oo 

Pon* = c*Y^ {~R^{Kt^'^)w*y c* ® nr. (11.15) 

n— 

Here, ("^ G Ranx/^^^o is the unique vector in the kernel of A^, normahzed as 

A#c# = o, IICII = 1, (CO-i- (11-16) 

(Recall that A* is either Aq or A^, so is either a null vector, of Aq or a null vector, (* , of Aq.) 
Letting C'^'^ equal C if # = * and letting it equal (* otherwise, the constant in (jll.lSp takes the 

form C* = (^C**(^ ^r, 

The overlap (C*iC) can be chosen strictly positive since the (^"^ are the Perron- Frobenius eigen- 
vectors of A"^ (i.e., their components can be chosen non-negative), and every component of ( is 
strictly positive (see below). The last relation in (|11.16|) is then achieved by scaling C* properly. 
The normahzation (fl*g,ng) = 1 together with (fTTTI) . pi.l2p (jll.lSp and (|11.16p imphes that 

C^C =l + o{g). (11-17) 

If the condition (|11.4p . r' ^ 5^'*'", does not hold then we have to apply the Feshbach map 
iteratively and use a corresponding perturbation theory for eigenvectors. We omit here this analysis 
and refer the reader to [HI El IS] for general references on such a RG perturbation theory and we will 
present elsewhere the RG perturbation theory in our specific case. 

In Section \W\ we have shown that NuUA = Cilp and consequently ( = flp and the vector 
f^o — C ® = ® is our Unperturbed state introduced in Section [S] Recall that Qp = ftp'' is 
the particle Gibbs state at temperature Pp. 

Expressions (|10.ip - (|10.2p for Ao imply the following relation among vectors (* corresponding to 
different particle temperatures 

C = ^(^^^{e^'"^'^ ® ]lp)CU,=o - (11-18) 

In view of (jll.lSp . it suffices to consider f3p — 0. 

In general there is no simple expression for the eigenvector C*. However, there are three cases 
where such an expression can be obtained. We expand C*|/3p=o in the basis ipj ® ifj as 

C\l3r,=0^^l3'Pj(^V3- (11.19) 
j 

Here cpj and Ej are defined after Eqn (|2.37p . Note that the normalization condition and Eqn (|2.37p 
imply that jj — ^fN . We have the following results: 
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(i) If /3i = /32 = /3 and (3p = 0, then 



(ii) If /?! is fixed and P2 — > 0, and /3p — 0, then 



7, = ^ + 0(/32). (11.21) 



(iii) If = dimTip = 2 and /3p = 0, then 

V2a(E) , V2 
a[liij + 1 a[-bj) + 1 

where 

Here, we use the notation E -.^ E2 — Ei, Pj{E) := -^-^ — - and 

g,{E):^ J [|G,(fc)i2|2 + |G,(fc)2iP] dfc. 

|fc| = |i5| 

Equation (|11.20p follows from (HHH), and (fTU:i| . The expressions (|11.22p ~ (|11.23l) come 

simply from solving a two-dimensional eigenvalue problem. Equation (jll.21[) follows from a straight- 
forward perturbation theory in ^2- See Appendix iP] for an outline of the proofs. 



12 Entropy production rate for rj 

In this section we prove Theorem l3.2l Recall that the stationary state 77 mentioned in this theorem 
is, in fact, given in (j9.6p . To analyze the entropy production, EP{r]), in this state rj we use expression 

EP{r,) = (/3i - /32)r/(0i), (12.1) 
which relates it to the heat flow, 77(01), in the state ry. Recall that 

01 = w[vi,Hr.i] = igi]{ai{ujGi) - al{ujGi)). (12.2) 

If /3i = P2, then ri{(f)i) — 0, c.f. [19]. We want to show here that ?7(0i) > if (3i > f32- Our proof 
is based on 

Theorem 12.1 Set (3 = max(/3i,/32) and let fl^ ^ (* with A*C* =0 (see U1.16\) ) and with 

the vector fig defined in \2.3T^ . Under the conditions of Theorem \3.1\ we have 

viM^g'v' + o{g')0{S(3), (12.3) 

where, recall, 5 [3 :— [3i — (32 and 

r\ = 2Re l?Ll, 7r(wi)L,.ii(io + i^)'^ e'^^^^l'^-Kivx)^^ . (12.4) 
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Moreover, we have the explicit expression 



j^Efae^-"--70 "-g^;^ , (12.5) 



where Eji ~ Ej — Ei, gji(E)'^ = J^j, (Pk\ {ipj,Gi{k)(pi) \'^S{Eji — lu). The numbers > are the 
components of the vector C*, see (|11.19p . normalized as in ill. 16]) . at j3p — 0. Observe that by 
mm . r]' = Ofor f3i=p2- 

The following result shows that r/ is strictly positive for small nonzero temperature differences. 

Theorem 12.2 If 5 [3 — j3i ~ (32 > is small and either dim Tip = 2 or the coupling functions \2.4^ 
satisfy G — Gi — G2, then the linear in 5(3 part of rj' is > Q. Moreover, in the latter case, this part 

IS 

6(3 Z,{(3,) ^ E% g^Ejk? 

2 ZJ(3i + (3J2) ^ e/^is. - e^^^^ ' ^ ' 

^ ^ j>k 

where Zp{(]) = tie^^^'' is the particle partition function. 

Remarks. 1. In the general case, if Gi is close to G2 one deduces strict positivity of 77' in the 
linear term in S(^ by a perturbation argument. 

2. For an expression of 77' in the case dim Tip = 2 similar to (|12.6[) . see the proof of Theorem 
112.21 given below. 

Proof of Theorem ll2.H To simplify the exposition we restrict ourselves to the case r' ^ 5^^", 
for a = ^Vt7§, and /i > 3/2. Pick po = g^^^". First we prove an estimate on 77(01) which is rougher 
than (|12.3p and then we explain how to obtain (|12.3p . Recall that 

7j(c(>i)^{n^,7r{({>,)ene). (12.7) 

Let ftf — ("^ ^flr and take a number M so large that {g'^^°'[mm{T' , g'^)]^^)^'^ — o{g). Substituting 
expansions (|11.7p and (|11.15p into the r.h.s. of this expression, using that (ff^^Po ^^^)^ — o{g), for 
a > i, and using that (pi is proportional to g, we find 



where 



7?((/.i) = 77o + 77i+Rem + o(g^), (12. 



7/0 (^!5,7^(0l)^^o>, (12.9) 

r/i ~g{ni[lgRo{LoeMc(,i)e+Tr{cbi)eRo{Loe)Ie)]^Q) , (12.10) 

M 1 



Rem = J2 E5'+'((-^o(^o5)/|)'=(-i?Qj(i^o^^'')VF*)'"i7S, 

m+n>l kd=0 

7:{(bM-MLoe)Iey{-RQ^^{K^o'^)Wrno). (12.11) 

Here we replaced on the r.h.s the factor C*G = 1 + 0(17) by 1 (see (jll.lSp and ()11.17p ). Using the 
pull-through procedure and elementary estimates of the resulting integrals we obtain that 

Rem==o(g2). (12.12) 

Since (f>i is linear in creation and annihilation operators, see (jl2.ip . we have 

770 = 0. (12.13) 
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It remains to compute 771. Using 

(Pp„ - l)I*n* = 0{po) (12.14) 

( [25j Lemma 5.3), removing the spectral deformation and using that =^ g[7r(wi), iL^i], we 

obtain 771 = g^r/' + o{g'^), where 

f]' = - (n*, [TT{vi)iLri{Lo + iQ)-^ ~ I{Ln + iO)~^iLriTT{vi)] I^q) ■ (12.15) 

Next, note that the contribution of the U2-part of / to rj' is zero since the resulting expression is hnear 
in creation and annihilation operators for the first and second reservoirs separately. The contribution 
of the 7r(?;i)-part of / is also zero by the symmetry of (|12.15p . Hence we have rj' = A — B, where 

A = {n*,n{vi)iLri{Lo+iO)-'7r'{-f'^\vi))no), 

B = {ni,n'{Y^\vi))iLriiLo+iOr\ivi)nn). 

Using that 7r'(7*/2(wi))l7o = je(I^L„-L)/2 

7r(wi)ilo, where L is given after (|4.19p . and the fact that 
Je~"^/^7r(ui)r2o = 7r(i'^)ilo = 7r(wi)rio, we transform 

A = {ni,T:{vi)iLri{Ln + iO)"^e"'^^«/V(ui)17o) . 
We use the relations 7r'(7*/2(t;i)) = Je(''^«~-^)/2^(wi)e~(^^o"-^)/V and 

JiLriiLo + iO)" V(ui)Oo = -iLriiLo + iOy^e^^^'^n{vi)no , 
see also (|5.ip . to find that 

B = -{n*, Je('^^°+^)/27r(t;i)iL,i(io + iO)-^e-l^^°/\{vi)no) . 
Finally, since {Ju, Jv) ~ {u,v), J^Iq = Qq, LpQg = and {(3Lo — L)Q^ = we obtain 

B = ~{%,TT{vi)iLriiLo + iO)-^e-'^^'>/^7r{vi)Qa) = -A. 

Since ij' ^ A- B this gives ^TM- 

Collecting estimates ([TZl)) . (|12.12p . (|12.13p and rji = g^rj' + o{g^) we find that 

viM = gS + o{g^) (12.16) 

where ry' is given by (|12.4p . This proves the rough version of (|12.3p - (|12.4p . 

Before we refine estimate (jl2.16p let us show (jl2.5p . We expand the vectors fl^ and fio in the 
basis (fjj ® fir , (pjj = ipj ® (fj , 

n n 
Slo = ^ (^jV^jj ^r, f^o = ^ Ij'Pjj (12.17) 

with aj > 0, J2j '^^ ~ ^ ^^"^ 7j — J2j '^jlj = 1- Plugging the expressions in (|12.17p into the r.h.s. 
of (|12.4p . using (j2.35p in order to express ■k{vi) in terms of creation and annihilation operators, 
and af^^ 

tt{vi) = a£i(^l + pi Gu) + arii^/fhGii) + h.c, 

where pi = {e^'^'^ — 1)^^) pulling through the annihilation operators to the right and using that 
7r(wi) (or Gie) acts only on the first (left) factor in ipjj = ipj (g) ipj, we obtain 

V' = -2^aj7,Ini /" + pi)u; (^^j,GliHp - + lu + iO)-^e-'^^"'--^^+'^'^/^Gi^,) 
-piw (^ifj,Gl{Hp -Ej-uj + iO)-^e-^^"''-'^'~'^'^'^Glipj^ ^d?k. 
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Inserting the partition of unity 1 = J^j \ iVjl inner products on the r.h.s. we obtain 

furthermore 

rj' = 2nJ2(^^l^ / u;\{ipi,Gifj)f{{l + pi)S{Eji+u;)-pi5{Eji-uj)}(fk. 

Interchanging the labels in the sum of the first term and noticing that in the resulting expression 
the integrals vanish unless Eji > 0, i.e. Ej > Ei, or j > i, we arrive at 

j>i 

where gji{E)^ J^3 d'^k \ {ipj,Gi{k)ipi) \'^6{Eji — uj). Observing that, due to (jll.lSp . 0^7^ = 
iV-^/^7.:|/3p=o = N~^^^-f,, we arrive at (fT231) . 

Since rf = 0{5P) (the energy flow vanishes if /3i = P2) estimate (I12.16P is ineffective if (5/3 is so 
small that (5/3 = o{g^). However, with a little bit more work p2.16p can be upgraded to estimate 
p2.3p . We sketch a proof of this estimate without going into much detail. We begin with some 
notation. 

Consider the self-adjoint Liouville operator for equal reservoir temperatures /32 = /3i = /3, 

L = Lo+gi, (12.18) 

where 

i = n{v) -n'{v) (12.19) 

with vr = 7r|^j=^2=/3 ^-^d similarly for tt'. Define = /2^" ' '^here L*^'^ = Lq + 7r(t;). 

Since the entropy production does not depend on /3p we set from now on /3p = Pi. The operator 
-ftr|/3p=/3i=/32=/3 = L is selfadjoint, (|5.1ip . and hence 51^ — il-g. From (|12.7p we obtain 

r/(0i) = ((r!^-%),^(,^i)el]e) + (%,^(0i)9(r!e-ile)) 

%,Tricj)i)eng). (12.20) 



We consider the last term first. Recall that ipi — g[v, iHri] — g[vi, iHri] and therefore Tr{(pi) — Tr{(j)i). 
It follows that _ _ _ _ _ _ 

(r!g,7r(,^i)eri(,) = (%,^((/)i)er!e) = (1^, i(</)i)r!). (12.21) 

Note that the r.h.s. of (|12.2ip describes the heat flow into reservoir rl for the equal temperature 
system. Since the heat flows vanish individually in the equal temperature case we have shown that 
the last term in (|12.20p vanishes. 

To estimate the first two terms on the r.h.s. we use as before expansions (|11.7p and pi.lSp 
for = ^e,^^ and similar expansions (obtained by setting = /3i = /3 in (|11.7p and (jll.lSp ) 
for rig. As a result we obtain an expression for (|12.20p of the type (|12.8p - (|12.1ip but with some 
of the powers in Rem, (|12.1ip . replaced by the differences, e.g. {—Ro{Loe)Iey ~ (— ^o(io0)^e)' or 
{Rq {Kq^^)WY'' — {RQ^{KQ^^)W)"-\f32=/3i=i3. These differences are estimated by using a telescopic 
expansion, e.g., 

(-i?o(Loe)/e)' - (-i?o(ioe)/e)' 
= J2^-MLoe)Iey-H~MLoeWe-IeK-MLoe)Iey-', (12.22) 
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and then estimating the first type of the differences in norm while for the second type we do first 
the pull-through and contraction procedure and then estimate the resulting integrals. As a result 
we have 



?y('/'i 



no 



where rji is given in (|12.10p and 

'72 = g i^o, 



-g {rio,TT{(t)i)0RQ{LQ^g) 

+o{g^)0{5p) 
V1+V2+ o{g^)0{5(3) 



l0Ro{LQ,g)T:{4>i)g + 7r(0i)ei?o(Lo,e)/ei 



(12.23) 



(12.24) 



Since the contribution of the ■i;2-component of Ig is zero we can omit the tilde (~) in (|12.24p . Thus 
the expression for 772 coincides up to the sign and the substitution fig — * VLq with the expression 
(IT^JU)) for 771, i.e. 



V1+V2 = -9 - ^Q, [l0Ro{Lo^g)Tr{(l)i)e + Tr{(l)i)gRo{Lo^g)Ig] rio) 



(12.25) 



We proceed with the r.h.s. of (|12.25p exactly as we did above with rji alone in equation p2.14p . and 
we arrive at 



rjicjji) = 2g2Re (q* - Q^, 7r(t;i)zL,i(Lo + i0)-ie-^^°/2^(«i)f^o) + o(<?2)0(<5/3). 
(|12.3p follows by noticing that 

Re (^no, TT{vi)iLri{La + iO)~^e~^^'>^^n{vi)Qa) =0. | 

Proof of Theorem [THU Consider first the case dim Hp = 2. Using p23)) . pi.22p . pi.23p 

we obtain in this case 



27r 



V 



a{E2i) + l 



E21 521 (£'21) 

gftSai _ 1 



Next we expand e^^^^^ — a(i?2i) around (5/3 = /3i — /32 = to verify that in the two-dimensional case 
the linear in (5/3 term of rj' is strictly positive. 

Now we consider the case Gi — G2- We want to control the components 7^, appearing in the 
expression (112. 5p for 77'. To this end we employ basic analytic perturbation theory (in (5/3) for the 
matrix family M{5(3) := Aq{SI3)\ 13^=0, where we consider /3i to be fixed. Write Tj instead of TjQ, see 
([TU31) . According to Proposition [TOTI we have M{SP) = Mq + S/SMi + 0{SP^), where 

Mo = -7(6-^1^-/2 ^]l){ri(/3i)-t-r2(/3i)}(e^i^-/2^]l), (12.26) 
Ml = -7(e-^i^''/2®]l){i[//p®]l,r2(/3i)]-(5^r2)(/5i)}(e'5i^''/2®]l). (12.27) 

Let = J2j>oi^f^y '^j ■ normalization (C*,C) — 1 (where C, — ilp{Pp) is the particle Gibbs 

state) implies that {Q,C) = 1 and (CjiC) = 0, for j > 1. Solving the zero-oder eigenvalue equation 
MoCo = gives 



(12.28) 
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The first-order eigenvalue equation reads MiQ + MqCi = 0, which implies 

[TiiPi) + r2{p,)]ief'^"-/' (g> m = dp\p,T2{P)ie^''''^^ ® ]l)Co- (12-29) 

We use here that {e'^^^p/'^ ^ is in the kernel of Hp. Let := Ylj^^'^^'^'^Vj ® Vj- Since 

V2W{P) = we have {dpV2){m{P) = -r2(/3)(90*)(/3), so 

9/3|ftr2(/3)(e'3^-/2®]l)Co = C(/3i)a^|;3,r2(/3)vE'(2/3i-/3) 

= C(/3i)r2(/3i)(Hp0]l)«'(/3i), (12.30) 

where C(/3i) = Zp{Pp) / Zp{(3i +/3p/2). The r.h.s. of (|12.30p is a vector in the orthogonal complement 
of kerr2(/3i) = C*(/3i). Using this fact and ^[TTSU^ we solve p2:29ll for Ci*: 

Ci* = C{Pi){e-'''"''^^<g>l)[ri{/3i)+T2{f3i)]-^T2{Pi){Hp(S>l)^{j3i) (12.31) 
+C'{e-^'"''/^ (g) (12.32) 

where the constant C is determined by the normalization condition {(^^ , =0. From expression 
(jl2.5p it is clear that the term (|12.32p does not contribute to the value of rf (this is the same as 
saying that rj' = for (5/3 = 0). 

Under the assumption Gi = G2 — G we have ri(/3i) = r2(/3i) and the r.h.s. of (|12.3ip simplifies 
to an easy expression, which, when used in (|12.5p . yields (|12.6p . | 

Proof of Theorem 13.21 If Gi = G2 or if the dimension of the particle system is 2 then we 
have rj' = rjSf] + 0{{5l3)'^) with ^ > independent of 5(3. This follows from Theorem 112.21 Hence 
for g and 5(3 both small, but independently of each other, we have, by (|12.3p . EP{r]p^i3^) > 0, which 
is the statement of Theorem l3.2l ■ 

Remark. A stronger statement, mentioned in the second paragraph after Theorem 13. 2[ can be 
proved as follows. By an abstract result of [H], EP{r]) > 0. Therefore, due to p2.ip . ry((/)i) > for 
Pi > 02- Hence, due to (|12.3p . for g sufficiently small (depending on 6 [3 in general), 

ri' >0 for Pi>(32. (12.33) 

Next, the are analytic in (3i and /?2 separately, away from f3i = 0, P2 = 0. To show this, we 
proceed as follows. From the explicit form of the level shift operator Aq, given in Proposition ! 10.11 and 
equation (jlO.Sp and similar expressions for diagonal elements, we know that AJ is analytic separately 
in (Ji and (32, everywhere except for /3i = 0, /32 = 0. We also know that for each /?i,/32 nonzero, Aq 
has a simple eigenvalue at zero (since Aq does). It follows from the Kato-Rellich theorem that we 
can find a zero- eigenvalue eigenvector CJ', which is analytic in Pi, (32- Next, we have to normalize 
that vector s.t. its overlap with Hp is unity. This yields (* = {Q,V,p)~^Q. Since {Q,Vlp) cannot 
vanish, we have that (* is analytic. Finally, the jj of Theorem 1 1 2 . 1 1 are the components of , with 
Pp = 0. From (|11.18p we see that analyticity of the components of (* in Pi,P2 is true regardless of 
what Pp is. Hence, the are analytic separately in P\,P2, for Pi 7^ 0, /32 7^ 0. 

Analyticity of the 7^ and expression (|12.5p show that 

77' is analytic separately in Pi and P2, for Pi 7^ 0, /32 ^ 0. (12.34) 

Equations (|11.2ip and (fT2?5l) imply that 

77' > if Pi is fixed and P2 is sufficiently small. (12.35) 

Relations (|12.33p - (|12.35p imply 'q{(l)i) > if Pi> P2 for almost all values of (/3i, /?2) e (0, 00) x 
(0,00), in the sense that for fixed Pi G (0,oo), ?/ can vanish only for finitely many values of P2 in 
any bounded subset of (0,oo). The same holds for Pi and P2 interchanged. 
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A Proof of existence of dynamics 

In this appendix we prove the existence of dynamics (|2.12[) . Recall the definition of the operator 
L'-'-) Lq + gTr{v) and the one-parameter group o-*(i?) := e'*^**'i3e~'*^"^' , B € tt^A)". 

Proposition A.l Assume the operators Vn A satisfy 12.13]) . Then for any state ijj normal w.r.to 
Uq the integrands on the r.h.s. of 112. 12]) are continuous functions, the series is absolutely convergent, 
the limit exists and equals 

V'*(A)=Tr(pa*(^(A))), (Al) 

where p is the positive trace- class operator defined by i'iA) — Tr((07r(A)). In particular, is 
independent of the approximating operators. 

Proof. Let w„ e ^ be an approximating sequence for the operator v satisfying (|2.13p . We define 
the selfadjoint operators Ln'^ :~ Lq + gn(vn) on the dense domain V^Lq). Let the one parameter 
group (T*^^ on 7r(^) be given by 

<^{7i)\-t^) — e -tJe 

Set crQ(7r(A)) :— Tr{aQ{A)) and let ip be an WQ-normal state on A, i.e. 

i^iA) = Tr(p7r(A)) (A.2) 

for some positive, trace-class operator p on ?i of trace 1. Then using the definition Vn = 7r(v„) we 
find 

^'{K- («„), • • • [a*^ («„), 4(A)] ...]) = Tr(p[a*" (U„), • • • [a*^ (K), a* (A)] •.•]). (A.3) 

Clearly the r.h.s. is continuous in ti, . . . ,tm and therefore the integrals in (|2.12p are well defined. 
Since the r.h.s. of (jA.3p is bounded by (c||Ki||)™||crQ(j4)||, the series on the r.h.s. of (|2.12p converges 
absolutely. In fact, using the Araki-Dyson series 

= E / • • • / dt„, K"- iAvn)), ■ ■ ■ 

m=0 

■■■[alHn{vn)),ali7T{A))]---l (A.4) 



one can easily see that the series in (|2.12l) is nothing but the Araki-Dyson expansion of the function 
Tr(pCT*^^^(7r(A))). Thus we have shown that the r.h.s. of (|2.12p is equal to lim„^oo Tr(p(T*j^^ (7r(A))). 

Now, Vn converges to V strongly on the dense set C :— TLp ® TLp ® Span{7r(A)rio}, where A 
ranges over all polynomials in creation and annihilation operators a*{f), j — 1, 2, with / G Lq. This 
follows from (|2.13p and the relation 

IKK - V)7T{A)no\\^ = c^o(A*« - V*){Vn - v)A). (A.5) 

Hence Ln'^ converges to L^^^ strongly on C. The set C is a core for both I'n^ and L^'^\ (This can 
be seen by using the GJN commutator theorem, [T3], Theorem 3.1, by taking Y = A + TV -f 1 for 
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the comparison operator in this theorem, and by noticing that C is a core for Y. The latter fact 
follows from [30], Corollary 2 to Nelson's analytic vector theorem X.39.) It follows from Theorem 
VIII.25 of that Li^^ converges to L'^' in the strong resolvent sense as n — > cx). In particular, 
e'*ii" e'*^'" strongly, so Tr(pcrf^J7r(A))) TT{pa\n{A))) which, in particular, shows (|JT|) . | 



B Positive Temperature Representation 
B.l Gluing 

In this appendix, we represent the Hilbert space Ti. in a form which is well suited for a definition of 
the translation transformation. This representation is due to |18) . 
Consider the Fock space 

T ■.= TiL^{X X {1,2})), X = RxS^ (B.1.1) 
and denote x = (u, a) X . The vacuum in J- is denoted by fir- The smeared-out creation operator 



a*{F), F e L'^iX X {1,2}) is given by 



'{F) = y2 [ F{x,a)a*{x,a) 
„, Jx 



IX 

and analogously for annihilation operators. The CCR read 

[a{x, a), a* {x' , a')] = 5a.a'5{x — x'). 
Following [18j . we introduce the unitary map 

U : [.F(L2(IR3)) ^ ^ [T{L'^{m?)) (g, TiL^iR^^))] ^ T{L^X x {1, 2})) (B.1.2) 

defined by 

U {[nrl ^rl] «) [flr2 «) flr2]) := (B.l. 3) 

and 

[a*(/i) «> 1 + 1 «> a*{gi)] (g) 1(g) 1 

+Ig)lg)[a*{f2)<»l + l<»a*ig2)])u-^ := a*{f®g), (B.1.4) 
where, for x — (u, a) e X, 

r.^ 1/ N jufaiua), u>0, m 1 ^^ 

[f®g]{u,(j,a):=i (B.1.5) 
\uga(—ua), u < V. 

This map is extended to the Hilbert space H = g) J- in the obvious way. We keep the same 
notation for its extension. 

The operators L^i g) lr2 + Iri 'S) ir2 and Nri g) lr2 + Iri ^ Xr2 are mapped under U to the (total) 
free field Liouvillian and number operator given by 



Lf = dT{u) — y I a* {x , a)ua{x , a) , 

a 

N = dr{l)^y2 a*{x,a)a{x,a). 
„, Jx 
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Moreover, the interaction takes the form 

UIU-^ ^a*{Fi) + a{F2), (B.1.6) 
where Fi,2 e L'^iX x {1, 2}, B{Hp (g) Hp)) are exphcitly given by (x = (u, cr) e X = M x S'^) 
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and 



J^2(«,a,a) = ^^— ^ (B.1.8) 

^ (Ga(MCT) ® Hp - e-/3="/2e-5/3.«/2 ^ Op (G^T* (m(t) )) , u > 

-x/^ (g;(-wc7) ® Hp - e-'3°"/2e-*/5°"/2 ]lp ® ap'^'^'^\'a^{-u<j))'^ , u < 

where (5/3q, — (3a — (3, SPp = (3p — (3, and /? = max(/3i, /?2)- 
Thus the operator K := UKU^^ can be written as 

K = Lo + gi 

where / — UIU^^ is given in (|B.1.6[) and Lq := ULqU^^ is of the form 

Lo = ip ® 1/ + Ip ® i/. 

B.2 Complex Deformation 

Now we express the complex deformation operators Ue introduced in Section [S] in the glued Hilbert 
space. For a function F £ L? {X x {1, 2}) and = ((5, r), x = {u, a) E X, define 

[ugF] iu,a,a) = ei''^''^^^F{jg{u),a,a), (B.2.1) 

where 

and sgn is the sign function, sgn(M) = 1 if w > 0, sgn(— u) = — sgn(u). Next, we lift the operator 
family ug from L'^{X x {1, 2}) to the operator family, Og, on (g) T{LP'{X x {1, 2})) in a standard 
way (cf. (|6.ip ). The family Ug is related to the family Ug introduced in SectionlHlas 

Ug = U-'UgU. 

The operator K :— UKU^^ becomes after spectral deformation 

kg := ilgKlJg^ = U,g + gig (B.2.3) 

where 

Zo,e = Lp + cosh (5 L/ + sinh (5 A/ + riV, (B.2. 4) 
A/ = dr(|u|) = >, / a* [x ^ a)\u\a{x , a) , 

ig = a*(Fi,e) + a(f^2,e) with Fj^g^ugF,. (B.2.5) 
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This spectral deformation can be translated to the original space Ti. as 

Kg := U^^kgU = Lo,9 + gle (B.2.6) 
where Lq^s :— U^^LQ^gll is given by (|6.14p and 

le = U-'ieU. (B.2.7) 

C The C*-algebra Ai 

Proposition C.l Let the algebras Ai and A be defined by \2.11^ with given in the proof 

below, and \2. 7\ ), respectively. Then Ai is strongly dense in A. 

Proof. Let = {/ a.e. continuous, ||/||^2 / \f{k)\'^{\k\^^ + l)d^k < oo}. Recall the real 

linear map 7^ defined in p.7p . Using that for 2: > 0, max(l/a;,l) < I'^l-a: < 4max(l/a::, 1), we 
obtain 

c||/||i. < WlpfWh = du -— \f{ua)\' < C\\f\\l. (C.l) 

so 7/3 is bounded and invertible (on R :— "jpLf C L^), and 7^^ is a (real linear) bounded map. We 
have the equivalence 

g e R ^ g{u, a) = -e''"/^g(-u, a) for a.e. u e M 

and y y da\g{u, a)]"^ < 00. (C.2) 

Let i?o := {g e i?, e^" g ^ for some 6 > 0} C i?. The set i?o is dense in R and (|C.2p implies that 

i?o = {e'3"/4/i| : e''"'/i € f^j. g^j^g 6 > and h{u, a) = -7^(-m, a)}. (C.3) 

Given g = e^^^'^h G i?0: define ft-e := * /i, the convolution in the variable u of h with the Gaussian 
Ge(u) := e~^G{u/e), where G{u) = 7r~^/^e~" and e > 0. /i^ is continuous (actually analytic), 
satisfies e**" /i^ € L^, and since G(:(-) is real valued and odd: 

h^{u,a) — —h^{—u,a). 

Therefore, <?e := e^^/^/i^ G i?o- 

Since ^ h in L^, we conclude that — > g in as e — > 0. Clearly, 
(?e extends to an entire function z 1— > gdz, cr). Define the set 

^anai _ {e/3«/4/j^| satisfies the conditions on r.h.s. of ((a3|) . e > 0}. 

jjanai jg gubget of i? that is dense in R. Since 7^^ is bounded, then — 7^^(i?^"''') is dense 

in L^. Since is dense in L^, we conclude that D'^'^'^^ is also dense in L^. 

Define Ai as in (|2.17[) with T)'^'^'^^ given above. Since T>^'^^^ is dense in L^, Ai is strongly dense 
in A (defined by (|2.7p V Next we have for real 9 

UeTT{Ap®Wifi)(^Wif2))Cl 
= TTpiAp)^W{i^0fi)e)®W{i^pf2)e)ne (C.4) 

where the map g ^ gg is defined by gg{u,a) := g{je{u),<j) with the function jg defined in (A2.2) 
and where we understand the Weyl operators on the r.h.s. as acting on the (glued) GNS space 
and Qg given in the same representation. Using that 9 {'lpfj)e are analytic for fj G 2?™'^' and 
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{^pfj)e S L^(R X S"^) as long as | tan(Im^)| < and expanding W{{jpfj)e) into the Taylor series 
and r^e into the Dayson-Araki one, one can show the r.h.s. of (jC.4p has an analytic continuation in 
9 into the neighbourhood of M given by 

{61 e C^||tan(Im^)| < and \lmT\<p-^]. 



D The vectors 

In this appendix, we outline the calculation of the vectors C* in the special cases mentioned in 
pi.l8p - (|11.23p . As mentioned after (jll.lSp . we may restrict our attention to f5p = 0. 

(i) (3i — — P- Proposition llO.il gives 

A* - [e-^^-/2 ^ ^ Jg/3i/,/2 ^ J ^ 

Therefore, [e^^p/'^ (g) ]l] ^* must be in the kernel of A|q + Ajq, which is spanned by fJ^. Thus 



The normalization is given by setting (fip'^'' °\C*) = 1- This yields the expression (|11.20p . 

(ii) Proposition llO.il and equation p0.3p imply that 

AS = zrio + 0(l) (D.l) 

where the operator 0(1) is bounded as (32 0. The matrix elements of Fio — rio(/32) in the 
basis {(fin (Xi (fin} are (see [5], Eqns (B21)-(B22)) 

(rio)m,n — ^m,n ^ ^ ^ ^rak (1 ^ran)^mni 

k—O.k^m 



where Smn is the Kronecker symbol, and where 

g/32|iJ^„|/2 r 

We expand Pio in (32'- 



rio = -^r' + o(i), (D.2) 



where 0(1) is bounded as /?2 ^ 0, and where 

(r')m,n = 27r(5„i_„ ^ |i?mfc| / C?0'|G2 (|i?mfe | , Cr)m_fe ^ 
— 2ti{1 — 5„in)\E„in\ / dcr|G2(|i?mn|, Cr) 

It is obvious that the vector with constant coordinates [1,1,..., 1]"^ is in the kernel of P', and, 
by the Perron- Frobenius theorem, that zero is a simple eigenvalue of P'. 

Hence, we see from (jD.l|) and (jD.2|) . by perturbation theory, that the vector in the kernel of AJ 

is of the form C* cx [1, 1, . . . , l]^+0(/32), as (32 0. The proper normalization {nf''^°\c*) = 1 
yields pT^ . 
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(iii) dim "Hp = 2. Let E = E2 — Ei > 0. We use Proposition 110.11 and formula ()10.5p to obtain the 
following expression for the level shift operator Aq: 



where 



Ao = 2TTiE^ 



a 

-b 



= E 



/ da\[G,iE,a)U 



S2 



da \[G,{E,a)U 



It is easily seen that flf'' 



[1, 1]^ is in the kernel of Aq, as it should be. The eigenvalues of 



Aq are thus zero and TrAg = 2TriE'^{a + b) ^ 0. The kernel of Aq is spanned by C[b/a,l] 
^* (X [b/a, 1]-^ = ^(pi 1^ (fii + (p2 (E) (p2- The normalization is given by setting {^lp^''~^\ (*) 
This yields the expressions (lll.22|) . (|11.23p . 



so 
= 1. 
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